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ABSTRACT 


Engineering  decisions  must  be  made  in  the  presence  of  uncertainties  which  arise  as  a 
consequence  of  imperfect  information  and  knowledge  and  inherent  randomness  in  many  design 
parameters.     It  is  on  account  of  these  uncertainties  and  potential  risks  arising  therefrom 
that  safety  margins  provided  by  load  and  resistance  factors  are  required  in  design.  Reli- 
ability methods  are  employed  in  this  study  to  facilitate  the  selection  of  criteria  for 
reinforced  concrete  design.     These  methods,  which  are  based  on  probability  theory,  provide 
a  logical  basis  for  determining  the  manner  in  which  uncertainties  in  resistance  and  loads 
affect  design  safety  and  how  their  effects  should  be  controlled  in  building  standards. 
Following  a  comprehensive  analysis  of  uncertainty  measures,  safety  indices  associated  with 
existing  reinforced  concrete  designs  are  evaluated.     Design  criteria  commensurate  with 
levels  of  uncertainty  and  required  reliability  are  then  presented.     Simplification  of  these 
leads  to  practical  reliability  based  criteria  which  retain  the  relatively  simple  charac- 
teristics of  existing  criteria  and  yet  have  a  well  established  and  documented  rationale. 

KEY  WORDS:     Buildings  (codes);  concrete  (reinforced);  design  (criteria);  loads;  probability 
theory;  reliability;  statistical  analysis;  structural  engineering. 
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NOTATION 

In  addition  to  standard  American  Concrete  Institute  notation,  the  following  symbols  are 
frequently  employed  in  this  report: 

A  =  influence  area,  equal  to  2A^  for  beams  and  4A^  for  columns; 

A^  =  tributary  loaded  area; 

B  =  bias  in  mathematical  model; 

b  ,  h  ,  d  ,  A      =  nominal  diminsions  for  width,   thickness,  effective  depth  to  reinforcement, 
n'     n'     n'     sn  ,       .   J,       .       ,  ^.  i 

and  reinforcing  bar  area,  respectively; 

E[X],  Var[X]  =  mean,  variance  of  random  variable  X; 

F      f    =  cumulative  distribution  function,  probability  density  function  for  variable  X; 

A  A 

k£/r  =  slenderness  ratio; 

D,L,W  =  Dead  load,  live  load,  wind  load,  respectively; 

L  ,       ,  L    =  code  -  specified  basic  live  load  instantaneous  live  load,  and  lifetime  maximum 
c      I  t 

live  load; 
p^  =  probability  of  failure; 

p..  =  relative  likelihood  of  load  ratio  combination  ii ; 

Q  =  maximum  variable  load  which  is  added  to  other  loads  taken  at  their  instantaneous  values 
when  combining  loads; 

R^,  Sj  =  resistance  variable  i,   load  j,  respectively; 

R*,  S*,  =  the  values  of  resistance  variable  i  and  load  i  which  are  substituted  in  the  limit 
11 

state  equation  for  safety  checking; 
R*  =  factored  ultimate  resistance; 
s  =  spacing  of  stirrups  used  to  resist  shear; 
T  =  time; 

U*  =  factored  ultimate  load 

X,  a    =  mean,  standard  deviation  of  variable  X; 
X 

=  direction  cosine  (see  Equation  2.8); 
3  =  safety  index; 

6  ,  A  ,  n    =  basic  variability,  prediction  error,  and  total  uncertainty,  respectively, 

A         X  X  . 

in  X; 

A",  A",  A"  =  tolerances  to  be  subtracted  from  b  ,  h  ,  or  d     for  safety  checking  when  using 
bhd  ^  ■   1   c     ^         r        ^  nn  n 

the  partial  factors  format; 

p,  p   ,  p     =  reinforcement  ratios  for  flexure,  shear  and  combined  bending  and  compression, 

^      ^      respectively;  p  =  A  /bd;  p  =A  /bs     ;  p     =  (A    +  A')/bh. 

s  V    V       St       g  s  s 
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<))',  y'  =  resistance  and  load  factors  (i  =  D,L,W)  applied  to  nominal  resistance  and  load; 

<j),        =  resistance  and  load  factors  applied  to  mean  resistance  and  load; 

Apostrophe  (prime)  denotes  a  nominal  value  specified  for  use  in  a  conventional  standard; 

e.g.,  nominal  compressive  strength  of  concrete,  f';  nominal  live  load 
L'  specified  by  ANSI  A58.1  -  1972; 

Asterisk  denotes  a  factored  value  of  a  variable  which  is  used  in  safety  checking. 
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1.  INTRODUCTION  AND  BACKGROUND 

The  design  of  safe,   functionally  reliable  and  economical  structures  is  of  continuing 
i|  concern  to  structural  engineers.     In  viev.-  of  existing  and  projected  shortcigrs  in  builHirp 

materials  and  continuing  increases  in  construction  costs,   those  costs  which  are  attributable 

i, 

I         to  excessive  safety  margins  can  and  should  be  reduced.   In  many  situations,  material  savings 
I  may  be  realized  without  sacrificing  structural  performance  and  safety. 
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Engineering  decisions  must  be  made  in  the  presence  of  uncertainties  which  are  invariably 
present  in  practice  [3,4]*.     These  arise  as  a  consequence  of  imperfect  information  and 
knowledge,  lack  of  experience,  and  inherent  randomness  in  many  design  parameters.     In  fact, 
it  is  on  account  of  these  uncertainties  and  the  potential  risks  arising  therefrom 
that  the  safety  margins  provided  by  load  and  resistance  factors,  and  the  like,  are  required 
in  design.     However,   the  safety  factors  provided  in  existing  codes  and  standards  have  been 
derived  primarily  from  experience,  tradition,  and  engineering  judgment  [51].     As  a  consequence, 
current  codes  may  differ  in  the  levels  of  safety  they  provide  for  different  materials  and 
modes  of  failure. 

In  recognition  of  these  factors,  the  current  trend  in  code  development  in  the  United 
States,  as  well  as  Europe,   is  toward  the  use  of  approaches  which  employ  probabilistic 
concepts  as  a  basis  for  developing  design  criteria.     This  represents  an  attempt  to  pro- 
portion members  and  connections  so  that  the  probability  of  their  failure  is  acceptably 
small.     In  modern  reliability  analysis,  probability  theory  furnishes  the  logical  framework 
for  incorporating  uncertainty  and  the  consequences  of  failure  which  are  reflected  in  the 
probability  of  failure.     The  safety  factors  obtained  are  functions  of  precisely  what  motivates 
and  necessitates  their  use  in  design,  and  are  commensurate  with  the  measure  of  uncertainty 
and  acceptable  failure  probability. 

An  example  of  reliability  based  criteria  is  given  by  the  tentative  load  and  resistance 
factor  design  (LRFD)  procedure  developed  by  Galambos  and  Ravindra  [28]   for  steel  structures. 
A  preliminary  study  of  such  criteria  for  reinforced  concrete  structure  was  conducted  by 
Ellingwood  and  Ang  [21].     In  Canada,  Allen  [1]  analyzed  the  reliability  of  steel  designs 
implied  by  using  the  load  factors  given  in  Section  4.1  of  the  National  Building  Code  of 
Canada.     Siu,  Parimi  and  Lind  [58]  carried  this  one  step  further,  in  determining  an  appro- 
priate load  and  resistance  factor  set  for  several  materials.     In  their  study,  and  in  the 
Canadian  Standard,  it  is  proposed  that  one  set  of  load  factors  be  used  irrespective  of  the 
material. 


^Numbers  in  brackets  denote  references  at  the  end  of  this  report. 
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In  Europe,  the  Comite  Europeen  du  Baton  (CEB)  has  also  proposed  a  material-independent 
loading  criterion  for  ultimate  and  serviceability  limit  states*  which  would  have  a  proba- 
bilistic basis  [12,26].     The  desirability  of  using  material-independent  loading  criteria 
appears  to  be  widely  agreed  upon.     This  represents  an  attempt  to  avoid  the  Inevitable 
conflicts  which  may  arise  when  different  material  usages  are  being  weighed  or  combined,  as 
in  composite  construction,  and  an  attempt  to  place  all  materials  on  an  equivalent  footing 
from  a  design  safety  point  of  view. 

It  should  be  pointed  out  that  the  Canadian  and  European  building  codes  prescribe  live 
and  wind  loads  that  may  differ,  in  some  cases  considerably,  from  those  presented  in  American 
standards  such  as  ANSI  A58. 1-1972  [9].     Therefore,  while  the  methodologies  employed  in 
those  studies  are  useful,  the  specific  criteria  and  numerical  values  cannot  be  assumed  to 
be  applicable  to  American  building  standards. 

In  the  United  States,  safety  criteria  have  traditionally  been  developed  separately 
for  different  materials  and  methods  of  construction  by  interested  groups  such  as  the 
American  Concrete  Institute  (ACI) ,  the  American  Institute  of  Steel  Construction  (AISC) , 
and  similar  organizations.     Needless  to  say,  this  has  led  to  a  situation  where  there  are 
many  differences  in  philosophy  and  inconsistencies  in  the  levels  of  safety  and  performance 
accorded  different  structures,  depending  on  the  materials  used  in  their  design.     A  stage 
Is  rapidly  being  reached  where  the  implementation  of  probabilistically  derived  design 
criteria  is  considered  both  feasible  and  desirable  in  the  United  States.     To  this  end, 
ANSI  A58.1  has  designated  a  subcommittee  to  investigate  the  possibility  of  developing  one 
set  of  load  factors  which  would  be  appropriate  for  use  with  different  materials.     In  order 
to  do  this  in  a  way  acceptable  to  current  material  groups,  appropriate  target  reliability 
levels  first  must  be  identified  for  different  construction  materials  and  design  criteria 
formats  must  be  decided  upon. 


*Limit  states  describe  various  conditions  where  the  structure  fails  in  some  way  to  achieve 
its  intended  purpose.     Ultimate  limit  states  relate  to  safety  under  maximum  loads,  while 
serviceability  limit  states  relate  to  loss  of  function  under  normal  loads. 
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The  objectives  of  the  present  study  are  actually  threefold.     The  first  is  to  provide 
a  convincing  demonstration  of  modern  reliability  analysis  techniques  in  selecting  practical 
design  criteria,  taking  advantage  of  recently  developed  [14,26,30]  methods  of  reliability 
analysis.     For  illustrative  purposes,  the  ACI  Standard  318-71  [10]  which  governs  much  of 
the  concrete  construction  in  the  U.S.  was  selected  as  a  framework  for  criteria  development. 
A  second  objective  is  to  determine  levels  of  safety  tacitly  implied  by  existing  reinforced 
concrete  design  standards.     This  then  contributes  to  the  technical  data  base  which  would 
be  required  prior  to  developing  material- independent  loading  criteria.     Although  the 
reliability  analysis  methods  used  in  this  study  differ  from  those  used  by  Galambos  in  his 
study  of  the  AISC  specifications,    [28,60]   the  numerical  results  obtained  in  these  two 
phases  may  be  considered  to  complement  his  work.     A  final  objective  is  to  examine  some  of 
the  relative  advantages  and  disadvantages  of  some  of  the  design  criteria  formats  most 
likely  to  be  proposed  for  use  in  first  generation  probability  based  standards. 

The  report  is  organized  to  show  how  a  large  aggregate  of  engineering  data  and  judgment 
can  be  systematically  combined,  using  probabilistic  techniques,  to  arrive  at  rational 
design  criteria  which  need  not  be  overly  complex.     Section  2  contains  a  discussion  of  the 
recently  developed  first-order  second  moment  reliability  analysis  method  which  serves  as 
the  basis  for  the  safety  checking  and  design  criteria  presented  in  later  sections.  The 
reader  familiar  with  the  mathematics  m.ay  proceed  directly  to  Section  3  without  signifi- 
cantly losing  continuity.     Section  3  contains  a  brief  discussion  of  how  the  measures  of 
uncertainty  can  be  estimated  and  a  summary  of  uncertainties  in  resistance  and  load  terms. 
The  analysis  of  Individual  uncertainty  measures  is  documented  in  detail  in  the  Appendix. 
The  preceding  tvzo  sections  lay  the  groundwork  for  Section  A,  in  which  reliabilities  associ- 
ated with  existing  reinforced  concrete  design  criteria  are  evaluated.     On  this  basis,  the 
development  of  reliability-based  design  procedures  are  discussed  in  Section  5;  practical 
load  and  resistance  factor  design  criteria  are  presented  in  Section  6. 


2.  RELIABILITY  ANALYSIS 

The  conceptual  framework  for  probability  based  design  is  provided  by  the  classical 
reliability  theory  described  by  Freudenthal  and  others  [27,4].     Tne  loads  and  resistances 
are  assumed  to  be  statistical  variables  and  the  necessary  statistical  information  is 
available. 
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A  mathematical  model  is  first  derived  which  relates  resistance  to  loads  by  describing 
the  limit  state  of  interest.     Express  the  overall  resistance  R  and  overall  load  effect  S  as 
functions  of  component  resistances  and  loads,  i.e. 

R  =        (R^,  R^.-.-R^)  (2.1a) 


S  =  gg  (S^,  S^.-.-S^)  (2,1b) 


where        and        are  (generally  random)  resistance  and  load  variables.     Since  a  structure  is 

safe  as  long  as  its  resistance  exceeds  the  effect  of  loads  placed  on  it,  safety  is  equivalent 
to  R  >  S,  or 

R  -  S  >  0  (2.2) 

Correspondingly,  the  limit  state  is  the  condition  R  =  S.  Then  safety  is  assured  by  assigning 
a  small  probability  p^  to  the  event  that  the  limit  state  will  be  reached,  i.e.. 


p^  -  P(R  <  S)  =  is)  fg  (s)  ds  (2.3) 

in  which  P(       )  is  the  probability  of  the  event  enclosed  in  parentheses,  and  F    and  f  are 

F  b 

cumulative  probability  distribution  function  and  density  function  for  R  and  S,  respectively 
If  the  probability  distributions  are  defined  by  their  means  and  variances  (or  coefficients 
of  variation,  the  square  root  of  variance  divided  by  mean).  Equation  2.3  may  be  expressed 
symbolically  as, 

Pf  =  Pf  (R,  s,  n^,  fig)  (2.4) 

where  R  and  fi    are,  respectively,  the  mean  and  coefficient  of  variation  of  R,  and 
R 

similarly  for  S.     The  c.o.v.  provide  a  convenient  dimensionless  measure  of  uncertainty. 


This  provides  a  theoretical  basis  for  the  evaluation  of  structural  safety.     It  is 

tacitly  assumed  that  all  uncertainties  in  design  are  contained  in  the  probability  laws  for 

R  and  S,  and  these  are  also  assumed  to  be  known.     However,  in  structural  safety  analyses 

these  probability  laws  are  seldom  known  precisely,  and  in  fact  only  the  first  and  second 

order  moments,  i.e.  mean  and  variance,  may  be  known  with  any  confidence.     While  this 

information  usually  is  sufficient  to  describe  the  central  behavior  of  the  probability 
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distribution,  it  may  be  grossly  inadequate  to  describe  the  behavior  out  at  the  tail  regions 
of  the  probability  distributions  which  are  significant  for  structural  reliability  analysis. 
The  problem  is  compounded  by  the  fact  that  the  probability  of  failure  computed  form  Eq . 
2.3  is  sensitive  to  the  choice  of  probability  distributions  for  R  and  S  at  probabilities 
considered  significant  for  structural  applications   [5].     Another  problem  is  that  while 
statistical  information  may  be  sufficient  to  define  R.  or  S .   in  Equation  2.1a  and  2.1b  the 
convolution  required  to  obtain  necessary  probability  laws  for  R  and  S  in  Equation  2.3  may 
be  impractical. 

The  difficulties  outlined  above  have  motivated  the  development  of  first  order  second 
moment  reliability  analysis  methods,  so  called  because  they  rely  on  a  knowledge  of  the  first 
and  second  order  statistical  properties  of  R  and  S  only  rather  than  the  probability 
distributions  as  a  whole.     Thus,  this  approach  is  generally  consistent  with  the  state  of 
knowledge  which  exists  in  structural  engineering  regarding  a  random  phenomenon.     The  concepts 
are  summarized  here,   there  being  several  recent  papers  on  the  topic  [3,14,26,30]. 

Let  the  limit  state  of  interest  be  given  by 


g  =  g^  (R^,  R^,...)  -  ggCS^,  S2,...)=  0 


(2.5) 


Failure  is  considered  to  occur  if  g  <  0.  For  purposes  of  design  or  safety  checking,  a 
set  of  points, 


R*  =  R.  exp  (-  a„     3        ),   i  =  1,  n  (2.6a) 

1  1  K ,  K . 

1  1 


=  S.    (1  +  3  fi^  ),  j  =  1,  m  (2.6b) 

>  ■'  3  3 


is  determined  such  that 


g(R^,  R2,...S^,  S^,...)  =  0  (2.7) 

An  exponential  transformation  is  used  on  resistance  terms  to  avoid  theoretical  difficulties 

*  *  * 

that  might  otherwise  arise  from  negative  R    values  [11,26].  Note  that  the  point  (R^,  R^ , . . . 

S^,  S2,...)  lies  on  the  failure  surface  defined  by  the  limit  state.     The  term  6  in 
Equations  2.6  is  the  safety  index  and  provides  a  measure  of  the  reliability  [30].     The  terms 
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a_     and  a„     appearing  in  Equations  2.6  are  the  direction  cosines  of  the  vector  normal  to 
R .  S . 

1  J 

the  failure  surface  g  =  0  along  which  3  is  measured  and  are  given  by, 

o  1/2 

=  a^/[Ea^^]  (2.8a) 

in  which 

a      =  a      exp  (-a      3  |f  ,  (2.8b) 

1  i  1  X        i '  * 

J  J       J  * 

where  o„    and  a„     are  standard  deviations  in        and  S.  and  the  derivatives  are  evaluated  at 
R .  S .  i  1 

1  2 

the  point  on  the  failure  surface  given  by  Equation  2.6.     The  reliability  analysis  requires 
that  Equations  2.6  -  2.8  be  solved  simultaneously  using  an  iterative  procedure  to  search  for 
the        which  minimize  the  distance  from  point  (R^,  R2»...i  S^,  S^,...)  to  the  limit  state 
equation . 

The  design  problem  is  to  find  an  appropriate  set  (R^^,  R^,...)  for  a  prescribed  safety 
index  and  load  set  (S^,  S 2- ••)  Conversely ,  safety  analysis  or  checking  requires  that  a 
value  of  3  be  found  corresponding  to  a  prescribed  set  of  resistance  and  load  parameters. 
In  either  event,  the  solution  is  a  function  of  the  reliability  3  and  the  measures 
of  uncertainty  Q. 


It  is  important  to  realize  that  if  the  probability  distributions  for  all  design 
variables  are  known,  then  the  safety  index  3  and  p^  in  Equations  2.3  and  2.4  are  uniquely 
related,  that  is,  for  any  prescribed  set  of  means  and  uncertainties  a  unique  3  corresponding 
to  p^  can  be  computed.     For  example,  if  the  limit  state  Equation  2.5  is  linear  in  the 
and  S.  variables  and  R^,        have  normal  probability  distributions,  then  3  =  *       (1  -  P^) , 
where  i""*"  is  the  percent  point  function  of  the  standard  normal  variate.     For  ease  of 
reference,  this  relation  is  shown  in  Table  2.1. 
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Table  2.1  -  Approximate  Relation  Between        and  B 

10"^  2.33 

10"^  3.09 

10""^  3.72 

10"^  4.26 


Even  when  such  distributions  and  probabilities  cannot  be  obtained,  however,  3  is  a  useful 
comparative  measure  of  reliability  and  can  serve  to  evaluate  the  relative  safety  of  various 
design  alternatives,  provided  that  the  uncertainties  are  handled  consistently.     This  will 
be  illustrated  in  detail  subsequently. 
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3.  MEASURES  OF  UNCERTAINTY  IN  RESISTANCE  AND  LOADS 


The  discussion  in  Section  2  has  shown  that  the  assessment  of  the  various  measures 
of  uncertainty  plays  a  central  role  in  reliability-based  design.     As  noted  previously, 
the  uncertainty  associated  with  a  particular  variable  is  measured  conveniently  in  a 
non-dimensional  form  by  its  coefficient  of  variation  abbreviated  c.o.v.     The  uncertainties 
used  in  the  reliability  analysis  should  include  all  imponderables  which  may  affect  design 
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safety,  and  would  include  not  only  inherent  statistical  variability  but  also  sources  due 
to  errors  in  prediction  and  in  imperfect  or  incomplete  information  [3,15].     The  key  test  in 
differentiating  between  these  two  sources  of  uncertainty  is  whether  the  acquisition  of 
additional  information  would  materially  reduce  their  estimated  value  [15].     If  the  variability 
is  intrinsic  to  the  problem,  additional  sampling  is  not  likely  to  reduce  its  estimate 
(although  the  confidence  interval  on  the  estimate  will  contract).     In  contrast,  uncertainties 
due  to  modeling  and  prediction  should  decrease  as  improved  models  and  additional  data  become 
available. 

In  this  section  general  methods  for  analyzing  uncertainties  are  discussed  and  results 
of  this  analysis  are  summarized  for  reinforced  concrete  member  capacities  and  structural 
loads.     A  detailed  examination  of  individual  measures  of  uncertainty  is  deferred  to  the 
Appendix  to  this  report. 

3.1  Methods  of  Analysis 

Let  X  denote  a  basic  resistance  or  load  design  variable.     The  true  statistical 
characteristics  of  X  should  be  employed  when  evaluating  p^  or  3    and  when  deriving  appro- 
priate factors  of  safety  for  design.     However,  the  random  characteristics  of  X  usually  are 
not  known  precisely  in  structural  engineering  problems  owing  to  insufficient  data  and 
imperfect  information.     Generally  available  instead  are  estimates  X  and  6^  of  the  mean 
and  c.o.v.  of  X.     These  are  usually  computed  from  data  gathered  under  ideal  or  carefully 

controlled  conditions.     Therefore,  while  6     reflects  basic  statistical  variability,  it 

X 

fails  to  encompass  all  sources  of  uncertainty  that  contribute  to  the  total  variability  in  X, 

9,  .     Portions  of  Q    may  arise  from  factors  which  are  inherently  random,  but  which  were 
X  X 

not  practical  or  feasible  to  include  in  the  experiment  or  its  analysis.     For  example, 

material  strengths  which  are  commonly  determined  for  design  purposes  from  standard  laboratory 

tests    may  differ  significantly  from  their  strengths  in  situ.     Load  idealizations  whose 

effect  cannot  be  evaluated  objectively  are  an  additional  source  of  uncertainty.  Moreover, 

despite  a  comprehensive  effort  to  identify  sources  and  magnitudes  of  design  uncertainty, 

there  will  invariably  be  some  that  are  neglected  or  unknown. 

These  various  omissions  constitute  prediction  and  modeling  errors.     If  the  uncertainty 

measure  attributed  to  these  factors  is  given  by  the  c.o.v.   A^^ ,   then  according  to  procedures 

described  in  detail  by  Ang  and  Cornell,    [3,4,15]  Q.      is  evaluated  as, 

X 
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(3.1) 


2        2  2 

Frequently,  A    can  be  broken  down  into  several  parts,   in  which  case,  A„  =          +  A.  +  .   .  . 

A  A  1  Z 

It  is  implicit  in  this  formulation  that  A    measures  primarily  the  uncertainty  in  predicting 

A 

the  true  mean  of  X  by  X.     Often,  A    must  be  estimated  on  the  basis  of  engineering  experience, 

A 

including  the  engineer's  feel  for  the  completeness  and  adequacy  of  his  assessment  of  the 
design  problem.     Thus,   the  capability  for  the  exercise  and  incorporation  of  engineering 
judgment  is  not  entirely  removed.     This  is  attractive  in  engineering  reliability  analyses, 
where  experience  is  often  as  valuable  an  indication  of  field  peformance  as  data  taken  under 
artificially  controlled  test  conditions  that  may  not  be  representative  of  actual  practice. 

When  data  are  available,  X  and  6„  can  be  estimated  from  the  data  samples  using  classical 

A 

Statistical  analysis  techniques.     However,   it  is  frequently  the  case,  at  least  in  structural 
reliability  analyses,  that  very  little  data  are  available  to  supplement  the  engineer's 
judgement.     Occasionally,  he  may  only  be  able  to  estimate  the  range  over  which  he  feels, 
from  past  experience,   the  data  should  lie.     In  this  case,   the  c.o.v.  may  be  estimated  from 
the  range  if  some  assumption  is  made  regarding  the  shape  of  the  probability  distribution  as 
well.     For  example,  if  it  is  assumed  that  values  between  x^  and        are  equally  likely 
(X  is  uniformly  distributed) ,  then 

~  i  ""2  "  ""l  (3.2a) 
'/3  ^2  ^  ^1 

If  it  is  assumed  instead  that  values  in  the  midrange  are  more  likely  than  those  near  the 

extremes  (X  has  a  bell-shaped  probability  density)  and  that  roughly  95  percent  of  the 
values  fall  within        and  x^,  then 

X     -  X 

6^i-^— (3.2b) 
X      2  x^  +  x^ 

Similar  techniques  may  be  employed  to  estimate  A  ,  provided  that  information  on  the  range  of 

A. 

the  means  is  available. 

The  resistance  and  load  terms  are  usually  functions  of  other  basic  random  variables  and 

in  order  to  apply  Equations  2.5  through  2.8,  their  means  and  variances  (or  c.o.v.)  must  be 

determined  from  the  known  statistics  of  the  component  basic  variables.     This  may  be 
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accomplished  using  linear  statistical  analysis,  which  involves  expanding  the  function 
in  a  Taylor  series  about  some  convenient  value,   truncating  the  higher  order  terms,  and 
taking  mathematical  expectations  of  the  linearized  result.     If  the  expansion  is  taken  about 
the  means  of  the  component  resistance  variables,   for  example,   this  leads  to  [4], 


E[R]  -        (R^,  R^,.    .    .)  (3.3a) 

1  1      1  J  J  111 

in  which  the  derivatives  are  evaluated  at  the  point  of  expansion  and  P^^  is  the  correlation 

coefficient  between        and  R_,  .     Usually,   the  basic  variables  are  selected  so  as  to  be 
mutually  uncorrelated ,   so  that  the  second  term  in  Equation  3.3b  is  zero. 

There  may  be  certain  instances  where  it  may  be  appropriate  to  expand  about  some  other 
point  than  the  means.     For  example,   the  reliability  analysis  described  in  Section  2  involves 
expanding  the  limit  state  function  about  a  point  on  the  failure  surface. 

In  addition  to  variability  in  the  basic  variables, the  mathematical  model  used  to 
describe  resistance  R  (or  load  S)  may  be  biased  and  may  contribute  to  the  overall  variability 
in  resistance  or  load.     This  is  because  analyses  of  structural  resistance  are  invariably 
based  on  assumptions  and  approximations  which  may  not  be  strictly  valid.     For  instance,  let 
R  be  the  actual  resistance  and  R    be  the  modeled,  or  predicted,   resistance.  Using  a  simple 
linear  model,  R  and  R^  are  related  by 


R  =  BR^  +  e  (3.4) 


in  which  B  is  model  bias  and  e  is  a  zero-mean  error  term.     The  mean  and  variance  of  R  are 

u 


given  by. 


R  =  B  R^  (3.5a) 


2—2     2  2 

aj,=  B    a^  +  a^  (3.5b) 
u 


2 

in  which        is  a  measure  of  the  observed  variability  in  R    on  the  condition  that  all  parameters 

14 


—  2 

in        are  known  and  where  R^  and  are  mean  and  variance  in  R  ,  which  may  be  found  using 

_  u 

Equation  3.3.     If  B  =  1,   the  model  is  unbiased.     As  with  Equation  3.1,   Equation  3.5  emphasizes 

the  contribution  of  two  separate  components  to  the  variance  in  the  actual  resistance  R.  If 
2—2  2 

=  R^  a  ,   conforming  to  the  commonly  observed  tendency  for  scatter  in  R  to  increase  as 
—  2 

R     increases,  B  and  a    may  be  estimated  as  the  mean  and  variance  of  the  ratio  R/R  ,  which 
u  ^  u' 

2 

is  a  common  way  of  comparing  test  data  with  theoretical  predictions.     More  precisely,   o  is 

mathematically  equivalent  to  the  conditional  variance  Var   [R/R   |r    =  x] .     Defining  ft    =  o  /R 

u    u  R  R 

and  n„    =  o„  /R  ,  we  have 
R  R  u 

u  u 

"r  "  "r^  ^  -"^"^^^  "^^-^^ 

u 

from  Equations  3.5a  and  3.5b.     The  second  term  in  Eq .   3.6  is  seen  to  be  equivalent  to  the 
A  -  term  in  Eq.   3.1,  discussed  by  Ang  and  Cornell   [3,4]  when  applied  to  the  strength 
prediction  equation  itself. 

3.2  Summary  of  Resistance  and  Load  Uncertainties 

An  analysis  of  individual  measures  of  uncertainty  in  the  capacity  of  reinforced  concrete 
members  and  in  structural  loads  was  performed  using  the  concepts  described  in  Section  3.1. 
This  analysis  is  described  in  detail  in  the  Appendix  to  this  report.     The  results  are  merely 
summarized  in  this  section,  so  as  to  better  maintain  continuity  in  the  overall  presentation. 
Along  with  the  c.o.v.  in  resistances  and  loads,  the  ratios  between  mean  and  nominal  values 
R/R',  S/S'  are  also  presented,  as  this  information  is  required  to  perform  the  calibration  in 
Section  4.     Finally,  statistics  used  in  other  recent  reliability  studies  are  also  given;  the 
similarity  of  the  results  obtained  may  be  noted,  despite  the  sometimes  different  methods  used 
in  their  derivation. 

Resistance  and  load  statistics  are  summarized  in  Tables  3.1  and  3.2,  respectively. 
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Table  3.1  -  Summary  of  Statistics  for  Resistance 


Source 

Flexure 

Shear 

Compression 
and  Bendine 

R/R' 

R/R' 

R/R' 

This  study 

1.12 

0.13-0.16 

1.18 

0.21-0.23 

0.97-1.13 

0.14-0.21 

Ellingwood,  Ang  [21] 

1.16 

0.16-0.17 

0.18-0.24 

0.14-0.22 

Siu,  et.  al  [58] 

1.14 

0.15 

1.10 

0.21 

1.14 

0.16 

Table  3.2  -  Summary  of  Statistics  for  Load** 


Source 

Dead 

Maximum 
Live 

Instantaneous 
Live 

Maximum 
Wind 

S/S' 

^S 

S/S' 

S/S' 

"s 

S/S' 

"s 

This  study 

1.0 

0.10 

Fig.  A. 6 

0.26 

Varies 

Varies 

1.2 

0.31 

Ellingwood,  Ang  [21] 

1.0 

0.13 

Varies 

0.33 

Varies 

Varies 

1.0 

0.37 

*Siu,  et.  al.  [58] 

1.0 

0.07 

0.7 

0.30 

0.80 

0.25 

Galambos  [28] 

1.0 

0.06 

Varies 

0.24 

Varies 

Varies 

1.0 

0.26 

*Allen  [1] 

1.0 

0.07 

0.7 

0.30 

Varies 

Varies 

0.8 

0.25 

*Siu,  et.  al  [58]  and  Allen  [1]  based  their  estimates  of  S/S'  on  the  Canadian  Standard  [43]. 

**The  fig  terms  include  a  component  to  account  for  the  analysis  used  to  transform  the  load 
to  a  load  effect  (moment,  shear,  etc.),  as  discussed  in  the  Appendix. 
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4.  CALIBRATION  WITH  EXISTING  DESIGNS 


Target  reliabilities  or  safety  indices,  along  with  the  measures  of  uncertainty  sum- 
marized in  Section  3,  provide  the  basis  for  reliability-based  criteria.     These  may  be 
established  either  by  prescribing  a  priori  values  which  are  consistent  with  overall  societal 
objectives  or  by  calibration  with  existing  building  codes.     The  calibration  approach 
appears  appropriate  in  the  short  term  because  it  references  the  reliability  analysis  to 
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existing  designs,  relying  on  the  consensus  that  provisions  in  current  design  codes  are 
generally  adequate  from  a  safety  standpoint,  and  thus  avoids  the  need  to  specify  an  acceptable 
risk  explicitly.     It  allows  the  usefulness  of  the  reliability  analysis  techniques  to  be 
demonstrated  in  a  familiar  context  and  enables  certain  inconsistencies  in  current  criteria 
to  be  eliminated. 


The  calibrations  in  this  section  determine  the  3  values  which  are  tacitly  associated 
with  design  criteria  in  Chapter  9  of  ACI  318-71   [10]   for  the  ultimate  limit  states.  These 
criteria  are  given  by, 

f  1.4D  +  1.7L';    [ACI  318-71  Eq.  9-1] 
<*>'  R'     >      f    0.75  (1.4D  +  1.7W'  +  1.7L');    [ACI  318-71  Eq.  9-2]  (4.1) 

^  0.9D  +  1.3W';    [ACI  318-71  Eq.  9-3] 
where  the  resistance  factor  '1''  depends  on  the  particular  ultimate  limit  state  (flexure, 
shear,  etc.)  under  consideration.     Other  loads  may  be  substituted  for  W.  Calibrations 
are  performed  on  an  individual  member  basis,  as  this  is  the  approach  normally  taken  in 
proportioning  members. 


Calibration  according  to  Equations  2.5  through  2.8  requires  that  the  mean  resistance 
obtained  using  existing  design  criteria  be  computed  for  various  load  combinations.  The 
dead  plus  maximum  live  load  case  is  considered  to  be  of  fundamental  importance.     It  is 
virtually  impossible  to  calibrate  against  all  possible  design  situations  or  load  combinations, 
since  not  only  would  the  number  be  prohibitive,  but  in  many  cases  the  information  available 
may  be  of  insufficient  quality  for  a  meaningful  calibration  to  result.     Instead,  it  appears 
appropriate  to  thoroughly  analyze  design  cases  which  have  demonstrated  satisfactory  past 
behavior  and  where  the  information  is  sufficient  that  the  S  values  obtained  from  calibration 
can  be  viewed  with  some  confidence.     The  gravity  load  case  is  the  most  obvious  of  these, 
and  indeed  controls  in  many  practical  design  situations.     Calibration  against  the  gravity 
plus  wind  load  combination  will  also  be  briefly  considered  for  comparison. 

4.1  Gravity  Loads 


Safety  indices  associated  with  structural  designs  which  conform  to  Equation  9-1  of  ACI 
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318-71  are  shown  in  Figures  4.1a,  b,  c  and  d  for  flexure,  shear  and  axial  compression 
combined  with  bending. 


The  live  load  specified  by  ANSI  A58. 1-1972  [9]  has  been  selected  for  these  calibrations. 

The  ANSI  design  live  load  L'  consists  of  a  basic  uniform  load  L     (e.g.,  L  =  50  psf  (2394 

c  c 

2  2 
N/m  )  in  offices)  which  may  be  reduced  when        is  less  than  100  psf  (4788  N/m  )  and  tributary 

2  2 

area  A^  exceeds  150  ft     (13.9m  ).  Thus, 


L'  =  L 

c 


1  -  min  [0,0008A  ,  0.6,  0.23  (1  +  D/L  )] 

t  c 


(4.2) 


This  relation  is  illustrated  in  Figure  A. 6  of  the  Appendix  for  the  case  where        =  50  psf 
2 

(2394  N/m  ).     If  L^/D  <  0.62,  the  maximum  live  load  reduction  of  60  percent  may  be  applied. 
This  would  be  the  case  for  many  reinforced  concrete  structures  in  which  the  dead  load  is  a 

major  portion  of  the  total  load.     However,  the  allowable  reduction  is  limited  where  L  /D 

c 

is  large.     An  examination  of  L^/D  for  typical  occupancies  in  reinforced  concrete  buildings 
including  offices,  corridors,  schools,  hotels,  and  retail  establishments  showed  that  L^/D 
may  vary  from  0.25  to  about  1.6,  with  a  substantial  percentage  being  less  than  0,8.  On 
the  other  hand,  the  live  load  of  interest  for  analyzing  the  reliability  is  the  maximum 
value  which  occurs  during  the  useable  life  of  the  structure.     Its  statistics  are  given  in 
the  Appendix;  although  these  are  based  on  an  analysis  of  office  live  loads,  the  limited 
data  on  c.o.v.  in  live  loads  for  other  occupancies  appears  similar.     The  mean  value  is 
compared  to  the  ANSI  live  load  in  Figure  A. 6.     The  ratio  between  the  nominal  live  load  and 
mean  live  load  is  a  function  of  the  tributary  loaded  area  and  L^/D. 

This  discussion  shows  that  3  for  existing  designs  will  depend  on  L^/D  and  A^ .  However, 

the  reduction  in  live  load  permitted  by  ANSI  A  58.1-1972  is  not  consistent,  in  a  statistical 

sense,  with  the  actual  reduction  in  the  mean  live  load.     This  is  reflected  in  the  dip  in 

2  2 

the  safety  index  (Figure  4.1)  which  occurs  e.g.,  at  A^  =  750  ft     (o9.7  m  )  when  L^/D  = 
0.5.     Similar  results  were  obtained  by  Galambos  and  Ravindra  [28]   in  their  calibration 
work  with  AISC  specif iciations  [60]  which  govern  steel  design.     In  contrast,  no  such 
irregularities  occur  in  calibrations  performed  on  the  Canadian  Building  Code  [1,43,58] 
because  the  live  load  reduction  formula  therein  has  a  form  similar  to  Equation  A. 13  which 
is  consistent  with  the  statistical  analysis  of  the  data. 
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When  L^/D  >  0.62  the  point  of  minimum  g  shifts  toward  smaller  tributary  areas.  This 

is  caused  by  the  limit  0.23(1  +  D/L  )  on  the  maximum  allowable  live  load  reduction,  which 

c 

shifts  the  point  of  minimum  L'/L  to  smaller  tributary  areas,  as  shown  in  Figure  A. 6.  While 
the  higher  live  load  and  its  large  attendant  variability  causes  the  reliability  to  be  less 
at  small        where  the  above  limit  is  not  effective,  at  large  areas  the  limit  on  load 

reduction  actually  causes  the  3  to  increase  as  L  /D  increases. 

c 

The  sensitivity  of  3  to  variations  in        is  also  shown  in  Figure  4.1.     As  a  coarse 

R 

rule  of  thumb,  an  increase  of  0.01  in  n    causes  a  decrease  of  0.10  in  3.     It  may  be  recalled 

K 

from  Section  3.2  and  the  Appendix  that  Q     in  flexural  resistance  decreases  from  0.16  as 

K 

the  depth  of  the  member  increases  because  the  variability  in  bar  placement  d  decreases  for 
deep  members. 

It  may  be  observed  from  Figure  4.1b  that  in  all  cases  the  3  value  for  shear  is  less 
than  for  flexure,  despite  the  desirability  for  beam  and  slab  failures  to  occur  in  the 
ductile  flexural  rather  than  the  relatively  brittle  shear  mode.     This  conclusion  has  also 
been  reached  in  other  studies  [21,58].     The  reason  for  this  should  be  clearly  understood. 
While  the  c.o.v.   in  flexural  resistance  is  of  the  order  0.13  -  0.16,  the  c.o.v.   in  shear 
capacity  is  of  the  order  0.21  -  0.23.     On  the  other  hand,  the  present  ACI  code  accounts 
for  this  increase  in  variability  and  the  undesirability  of  a  shear  failure  by  subjectively 
reducing  the  resistance  factor  from  0.9  for  flexure  to  0.85  for  shear.     This  small  reduction 
is  inadequate  to  account  for  the  increased  variability  and  desired  reliability  in  shear. 
This  is  an  excellent  example  of  the  type  of  inconsistency  which  may  be  common  in  current 
building  codes  because  of  their  failure  to  consider  uncertainties  explicitly. 

Figure  4.1c  shows  that  the  safety  for  tied  columns  is  greater  than  for  beams,  which 
is  in  accordance  with  ACI  objectives;   the  corresponding  difference  in  3  implies  over  an 
order  of  magnitude  difference  in  risk.     As  shown  by  comparing  Figures  4.1c  and  d,  the 
increase  in  ()>'  to  0.  75  for  columns  with  spiral  reinforcement  results  in  a  decrease  in  3 
of  roughly  0.3  from  that  for  tied  columns  at  all  tributary  areas.     The  spirally  reinforced 
columns  behave  in  a  more  ductile  fashion  than  tied  columns  and  it  is  primarily  on  account 
of  this  ductility  that      is  increased  to  0.75.     Thus,  the  difference  in  3  is  apparently 
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what  additional  ductility  is  currently  worth  to  ACI,  in  terms  of  reliability.     The  ACI 
code  also  currently  permits  <t> '  for  axial  compression  combined  with  bending  to  be  increased 
toward  the  cj) '  for  bending  when  the  axial  load  is  small.     This  increase  in  <}> '  appears 
justifiable  primarily  on  the  grounds  of  the  associated  increase  in  ductility  (see  Appendix). 

For  purposes  of  comparing  reinforced  concrete  and  steel  designs,  the  same  analysis 

was  performed  for  steel  beams  designed  according  to  Section  2  of  the  AISC  specifications 

[60].     Therein,        =  1.0  while        =  Yj^  =  1.7.     According  to  Galambos  and  Ravindra  [28], 

R'/R  =  0.93  ^nd  Q    =  0.13.     The  corresponding  safety  indices  are  shown  in  Figure  4.2  for 
R 

various  L  /D,  which  tend  to  be  higher  for  steel  construction  than  for  reinforced  concrete, 
c 

It  appears  from  Figure  4.2  that  an  appropriate  B  would  be  about  3.25  for  the  design  of 
steel  beams.     Galambos  used  a  lognormal  type  of  reliability  analysis  suggested  by  several 
researchers  [3,4,52]  which  differs  from  the  one  described  in  Section  2.     A  portion  of  his 
results  are  also  shown  in  Figure  4.2  for  comparison. 

Several  conclusions  may  be  drawn  from  this  analysis  of  concrete  and  steel  flexural 
designs.     First  of  all,  a  comparison  of  Figures  4.1a  and  4.2  shows  that  current  gravity  load 
design  criteria  for  reinforced  concrete  and  steel  beams  appear  to  result  in  similar  levels 
of  reliability.     Secondly,  there  does  not  appear  to  be  a  practically  significant  difference, 
at  least  in  this  case,  in  the  results  obtained  using  the  two  different  reliability  analyses, 
as  may  be  seen  in  Figure  4.2    Finally,  the  use  of  one  overall  factor  of  safety,  such  as 
the  1.7  prescribed  by  AISC,  Section  2  (or,  for  that  matter,  as  implied  by  any  working 
stress  provision  where  the  factor  of  safety  is,  effectively,  the  reciprocal  of  the  ratio 
of  allowable  to  ultimate  stress)  leads  to  more  scatter  in  the  resulting  value  of  3  than 
the  use  of  the  split  load  factors  followed  by  ACI  318-71.     This  may  be  seen  quite  clearly 
in  Table  4.1. 


Table  4.1  -  Comparison  of  3  for  Steel  and  Concrete 
Designs  when  A    =  500  ft^(46.5  m^) 


L  /D 
c 

0.  25 

0.50 

0.75 

1.0 

1.5 

Steel  (fi„  =  0.13) 

3.731 

3.551 

3.357 

3.182 

2.984 

Concrete (Q„  =  0. 16) 

3.028 

3.004 

2.946 

2.881 

2.839 
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O  L /D  =  I.O  (From  Ref.  28,  FlgureS) 


/ 


500  1000  1500  2000 


Figure  4.2  -  Safety  Index  for  Steel  Beams  Conforming  to 

AISC  Specifications,  Section  2  (Gravity  Loads) 
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As  the  L^/D  increases,  the  larger  variability  in  L  becomes  increasingly  more  important  for 
reliability  and  must  be  accounted  for  by  a  larger  load  factor.     It  may  be  inferred  that 
the  split  load  factor  approach  is  more  desirable  than  the  single  load  factor  approach  in 
situations  where  one  overall  load  criteria  (independent  of  load  ratio)  must  be  given. 

4.2  Gravity  Loads  Combined  with  Wind 

Calibration  performed  on  structural  designs  conforming  to  Equation  9-2  and  resistance 

factors  of  Section  9.2  of  ACI  318-71  are  shovm  in  Figures  4.3a  and  4.3b  for  flexure  and 

thrust  in  tied  columns.     6  is  shown  as  a  function  of        and  is  computed  for  both  load 

combinations  in  Equations  A. 17a  and  A. 17b.     It  may  be  seen  that  Equation  A. 17b  represents 

the  critical  combination  for  the  case  illustrated  (L  /D=0.5).     3  is  virtually  constant 

c 

when  A^  is  large  because  the  ratio  L'/L^  is  constant  at  large  areas  for  short  term  live 

load  and  the  contribution  of  the  large  variability  Q     is  dampened  by  the  typically  low 

^I 

values  of  L^/D.     (As  seen  from  Figure  4.1,  this  is  not  the  case  for  maximum  live  load). 

An  example  of  safety  indices  achieved  by  existing  column  designs  is  given  in  Table 

2  2 

4.2.     The  column  supports  a  tributary  loaded  area  of  750  ft.     (69.7ra  ).     This  corresponds 

to  the  "worst"  case  of  Figure  4.2  in  that  the  ratio  L'/L^  is  smallest  at  this  area.  The 

slight  tendency  for  the  reliability  to  increase  with  increasing  L^/D  may  be  noted,  as  well 

as  its  relative  constancy  with  W/D  at  moderate  values  of  L  /D. 

c 

Table  4.2  -  Safety  Index  for  Existing  Tied  Column  Designs 


A^  =  750  ft^     (69. 7m^) 


L  /D 
c 

W  /D 

0 

25 

0.  50 

0. 

75 

1 

0 

1 

5 

0.5 

2 

680 

2.825 

3. 

084 

3 

417 

3 

907 

1.0 

2 

727 

2.839 

3. 

041 

3 

309 

3 

741 

2.0 

2 

698 

2.773 

2. 

911 

3 

100 

3. 

431 

3.0 

2 

667 

2.  724 

2. 

827 

2, 

971 

3. 

233 

=  0.7; 

Yj^  =  1.05; 

=  ^W  = 

1.275 
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4.3a  -  Flexure  (L  /D  =  0.5) 
c 
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4.3b  -  Compression  and  Bending  in  Tied  Columns  (L  /D  =  0  5) 
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Figure  4.3  -  Safety  Index  for  Designs  Conforming  to 

ACI  Standard  318-71  (Gravity  plus  Wind  Loads) 
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For  consistent  reliability,   it  would  be  desirable  for  the  6  values  to  be  the  same  for 
both  (D  +  L)  and  (D  +  L  +  W)  loading  combinations.     However,  comparison  of  Figures  4.1  and 
4.3  indicates  that  the  probability  of  failure  under  the  (D  +  L  +  W)  combination  is  ap- 
proximately one  order  of  magnitude  higher  than  for  (D  +  L) .     There  are  a  number  of 
reasons  why  the  actual  reliability  under  wind  load  may  exceed  the  apparent  values  in  the 
analysis  above.     Among  these  are  the  wind  directionality  effect  alluded  to  earlier,  the 
fact  that  the  extreme  wind  may  not  act  in  the  direction  most  unfavorable  to  structural 
response.     In  addition,  elements  which  are  not  designed  as  load  bearing  may,  nevertheless, 
carry  a  certain  portion  of  load  and  thus  reduce  the  wind  load  resisted  by  the  structural 
system.     Lateral  drift  limitations  for  tall  buildings  may  also  require  that  the  lateral 
load  carrying  system  be  stiffened  above  that  required  by  an  ultimate  strength  analysis  alone. 
There  remains,  however,  the  very  real  possibility  that  current  reinforced  concrete  designs 
may  not  have  as  high  an  implied  margin  of  safety  for  combined  loads  as  for  gravity  loads 
alone. 

When  lateral  and  gravity  loads  counteract  one  another,  as  e.g.,  the  case  for  columns 

on  the  windward  side  of  a  building,  designs  must  conform  to  Equation  9  -  3  of  ACT  318-71. 

2 

Corresponding  3  values  obtained  for  columns  with  tributary  areas  exceeding  750  ft 
2 

(69.7m  )  are  shown  in  Table  4.3. 

Table  4.3  -  Safety  Index  for  Tied  Columns  with  Counteracting  Loads 


A    =  750  ft  ^       (69. 7m^) 
t 


L  /D 
c 

W'/D 

0, 

25 

0.  50 

0.75 

1 

0 

1. 

5 

1.0 

2, 

680 

2.846 

2.993 

3 

126 

3. 

359 

2.0 

2 

754 

2.840 

2.920 

2 

995 

3. 

134 

3.0 

2 

727 

2.784 

2.839 

2 

891 

2. 

990 

=  1.3 
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Safety  indices  for  smaller  W /D  are  not  given,  as  one  of  the  two  preceding  load  combinations 
would  normally  govern  design  in  that  case.     Neglecting  the  beneficial  contribution  of  live 
load  when  gravity  and  wind  loads  effects  counteract  one  another  is  obviously  conservative, 
as  shown  by  the  increase  in  3  with  L^/D.     Comparison  of  the  results  in  Tables  4.2  and  4.3 
shows  that  Equations  9-2  and  9-3  of  ACI  318-71  are  basically  risk-consistent  for  moderate 

values  of  L  /D. 

c 

4.3  Target  Reliabilities  for  Design 


Target  safety  indices  for  reinforced  concrete  design  were  selected  considering  all  of 
the  calibrations  performed  in  the  previous  sections.     In  the  study  by  Siu,  et.  al.    [58]  on 
the  Canadian  code,  a  weighted  average  of  all  the  3  -  values  was  taken,  weighted  by  the 
probabilities  of  occurrence  of  the  various  mean  load  ratios  (L/D,  W/D) .     However,  the  live 
load  reduction  employed  in  the  Canadian  Standard  is  basically  consistent  with  statistical 
analysis  while  the  ANSI  reduction  formula  is  not.     As  a  consequence,  Siu's  values  of  B  do 
not  depend  on  the  tributary  or  influence  area  while  the  ones  in  this  study  do.  Thus, 
appropriate  target  3  values  were  necessarily  determined  more  subjectively  in  this  study. 

For  flexure,  the  value  of  3  selected  was  3.0,  for  tied  columns,  3.5,  and  for  spirally 
reinforced  columns,   3.2.     These  values  will  tend  to  provide  slightly  more  safety  than 
existing  designs  at  large  tributary  areas  and  slightly  less  at  small  areas.  Considering 
consequences  of  failure,  additional  safety  would  probably  be  desirable  at  large  areas.  The 
difference  in  3  implies  that  columns  will  be  designed  with  a  probability  of  failure  which 
is  at  about  one  order  of  magnitude  less  than  for  beams.     While  this  may  not  be  deemed  to  be 
a  sufficient  additional  margin  of  safety  in  view  of  the  relative  consequences  to  the  structure 
of  the  beam  and  column  failure,   it  appears  justifiable  on  the  basis  of  current  practice. 
This  points  out  one  disadvantage  in  relying  too  heavily  on  calibration  to  set  design  safety 
standards.     On  the  other  hand,  an  increase  in  3  to  4  or  more  for  columns  would  result  in 
designs  which  are  more  conservative  than  currently  acceptable  for  atl  A^  and  would  probably 
be  rejected  by  the  profession  on  economic  grounds. 
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The  value  of  3  for  shear  presents  a  problem  in  that  current  values  obviously  are  too 
low  and  are  incompatible  with  the  desired  objective  of  having  members  fail  in  flexure 
rather  than  shear.     In  this  study,   6  was  increased  to  3.0  for  shear,  this  being  the  minimum 
adjustment  necessary  to  make  the  reliability  in  shear  no  less  than  that  in  flexure.  While 
it  may  well  be  desirable  that  3  for  shear  exceed  3  for  flexure,   further  discussion  of  this 
point  by  standards  committees  appears  advisable  because  of  the  impact  that  such  a  recom- 
mendation would  have  on  current  design.     It  should  also  be  noted  that  3  for  columns  ranges 
from  3.2  to  3.5.     In  view  of  the  relative  consequences  of  beam  and  column  failure,  it 
certainly  does  not  appear  necessary  to  increase  3  for  shear  above  3.2. 


In  contrast,  Galambos  [28]  prescribed  3  =  3.0  for  all  structural  members  in  his  develop- 
ment of  load  and  resistance  factor  criteria  for  steel  design.     While  simplifying  criteria 
development,   the  relative  consequences  of  failure  in  beams  and  columns,  the  latter  of  which 
is  typically  more  severe,  thus  is  not  considered.     It  should  be  observed  from  the 
above  discussion  that  the  use  of  calibration  as  the  sole  basis  for  developing  reliability 
based  design  criteria  implies  that  the  current  code  provisions  are  internally  consistent. 
Thus,  inconsistencies  of  the  type  found  in  shear  would  tend  to  carry  over  into  the  reliability 
based  code  if  the  results  of  the  calibration  are  not  treated  judiciously. 
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5.  DEVELOPMENT  OF  RELIABILITY  BASED  DESIGN  CRITERIA 

Specific  design  criteria  are  developed  on  the  basis  of  target  reliabilities  furnished 
by  the  calibrations  discussed  in  Section  4.     Although  many  different  formats  are  possible, 
it  is  felt  that  simplicity,  consistency  or  continuity  with  existing  formats,  and  ease  of 
enf orcibility  initially  should  guide  the  selection  of  these  criteria.     Such  practical 
considerations  are  necessary  in  order  for  the  approach  to  be  accepted  by  the  profession 
at  large. 
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Three  different  levels  of  sophistication  for  reliability-based  design  can  be  identified. 
At  all  three  levels,  statistical  variabilities  in  resistance  and  loads  are  taken  into 
account,  however.     At  the  highest  level,  termed  Level  III,  the  probability  distributions 
of  all  variables  and  necessary  convolutions  must  be  known  or  computed,  and  the  design  is 
based  on  an  acceptable  small  probability  of  failure  previously  agreed  upon  by  the  profes- 
sion.    As  noted  before,  this  information  generally  is  not  available  at  present.     Level  II 
methods  employ  safety  checks  at  discrete  points  of  the  limit  state  equation  (e.g.,  at 
selected  values  of  L/D,  W/D,  etc.),  the  basic  design  variables  in  the  limit  state  equation 
being  specified  in  advance.     Probability  distributions  are  not  required  for  Level  II 
methods  in  which  reliability  is  measured  by  the  safety  index  3.     Level  I  methods  involve 
the  selection  of  one  set  of  load  and  resistance  factors  to  be  applied  to  all  designs, 
regardless  of  L/D,  W/D,  etc.;  current  design  criteria  are  of  the  Level  I  type.     Levels  I 
and  II  can  be  made  equivalent  if  the  load  and  resistance  factors  are  allowed  to  vary  in 
accordance  with  the  load  ratios. 

In  this  section,  Level  II  design  procedures  for  flexure,  shear  and  thrust  are 
presented.     The  criteria  in  this  study  are  presented  in  terms  of  mean  resistance  and  load 
values;   to  convert  to  the  nominal  resistance  and  load  values  that  designers  may  be  more 
familiar  with,   simply  multiply  the  factors  by  the  ratio  of  the  mean  to  nominal  value. 
(The  nominal  value  is  the  code-specified  minimum  value  used  in  design;  for  example,  for 
Grade  40  reinforcement,  the  nominal  yield  stress  is  40  ksi,  while  the  mean  may  be  about 
47  ksi,  depending  on  the  bar  diameter.) 

Specific  design  criteria  generally  take  the  form  of  an  equation  describing  some 
limit  state  in  which  certain  "resistance"  and  "load"  factors  are  appended  to  terms  on  the 
resistance  and  load  sides  of  the  equation,  respectively.     The  safety  check  for  ultimate 
limit  states  requires  that  the  factored  ultimate  resistance  R*  exceed  the  factored  load 
U*,  or 

R*  >  U*  (5.1) 
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There  are  at  least  two  ways  of  specifying  the  factored  resistance  R* ,     In  the  first, 
followed  by  the  current  ACI  provisions,  one  overall  resistance  factor  (RF)  4)  is  appended 
to  the  resistance  side  of  the  design  equation  for  each  limit  state.     Using  this  approach, 
the  factored  ultimate  resistance,  R*  for  a  particular  limit  state,  thus  becomes, 

R*  =  (1)  R  (5.2) 

Here  the  basic  resistance  variable  is  R  itself,  and  the  resistance  function  g     in  Equation 

R 

2.5  is  simply  R.     This  is  also  the  format  selected  by  Galambos   [28],   for  steel  design. 
It  has  the  virtues  of  simplicity  and  familiarity  to  American  designers,  as  well  as  an 
ability  to  incorporate  the  consequences  of  failure  in  a  particular  limit  state  in  the 
resistance  factor  for  that  limit  state. 

A  second  approach  employs  partial  factors  (PF)  which  are  applied  to  the  means  of  the 

individual  variables  in  the  resistance  equation.     For  example,  let  f*  =  *     f  ,  f*  = 

y       y    y  c 

d)  f  ,  h*  =  h        -  A."  and  d*  =  d        -  A",  where  A,"  and  A"  are  small  dimensions  (on  the 
c  c  n  h  n  d  h  d 

order  of  1/4  -  1/2  inch  (6.4  -  12.7  mm))  subtracted  from  nominal  dimensions  for  purposes 
of  safety  checking.     The  factored  ultimate  resistance  is  then, 


R*  =  R  (f*,  f*,  h*,  d*)  (5.3) 

y  c 


In  this  approach,   the  partial  factors  4)   ,  <}>   ,  etc.  are  specified  only  once  in  the 

y  c 

design  standard  and  the  same  factors  would  be  used  for  all  ultimate  limit  states.  This 
format  was  discussed  by  ACI  Committee  318  over  a  decade  ago  [48] ,  and  has  become  popular 
in  Europe  [12],     The  application  of  safety  factors  directly  to  the  parameters  which  are 
the  source  of  variability  is  an  attractive  feature.     Despite  certain  conceptual  advantages, 
this  approach  was  not  acceptable  tc  many  American  engineers  and  was  rejected  by  ACI  member- 
ship when  it  was  originally  proposed,  as  there  was  a  feeling  that  it  might  penalize  concrete 
construction  in  relation  to  steel.     This  point  will  be  discussed  further  after  specific 
reliability-based  partial  factors  are  derived. 


Similarly,  the  load  side  of  the  equation  may  be  treated  several  ways.     The  simplest 
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is  to  apply  one  load  factor  to  the  sum  of  whatever  loads  the  structure's  limit  states 
must  be  checked  against;  the  factored  ultimate  load  would  then  be. 


U*  =  Y  (D  +  L  +  •  •  •)  (5.4) 

This  is  equivalent  to  AISC  Section  2  [60] .     This  format  recognizes  that  it  is  the  sum  of 
the  loads  acting  on  the  structure  at  any  time  instant  which  is  important.    However,  there 
has  been  a  tendency  to  attach  separate  load  factors  to  each  load.    Among  the  reasons  for 
this  are  that  the  variability  in  permanent  and  transient  loads  are  often  of  different 
orders  of  magnitude  and  that  individual  design  loads,  such  as  live  load,  wind  load,  and 
earthquake  load,  are  usually  determined  by  different  technical  committees.     The  factored 
ultimate  load  in  this  format  is 

U*  =  Y^D  +  Y.L  +  Y„W  (5.5) 
D  L  W 

or  some  variant  of  this  equation. 

In  this  study,  criteria  of  the  form 

(1)R  >  YqD  +  Y^L  +  Y^JW  (5.6) 

will  be  emphasized  because  it  is  easiest  to  tie  these  to  existing  ACI  criteria  and  the 
criteria  proposed  by  Galambos  [28]  for  steel  structures.    However,  the  other  criteria 
will  also  be  examined.     It  is  not  the  intent  here  to  advocate  any  particular  one  of  these 
alternatives;  indeed,  each  has  certain  advantages  and  disadvantages.     Regardless  of  which 
format  is  selected,  however,  the  safety  factors  determined  by  the  reliability  methods  are 
correctly  apportioned  to  the  resistance  and  load  sides  of  the  design  equations. 

It  is  assumed  in  this  study  that  the  total  load  effect  is  a  linear  combination  of 
the  individual  loads.     The  limit  state  equation  may  then  be  expressed  in  terms  of  the 
mean  load  ratios  by  dividing  Equation  2.7  by  D;   thus,  the  importance  of  L/D,  W/D,  etc.  in 

Sections  5  and  6.     The  load  factors  in  Equations  5.4  through  5.6  should,  in  general,  be 
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applied  to  the  load  prior  to  performing  the  analysis  which  transforms  the  load  to  a  load 
effect.     This  is  also  required  by  ACI  318-71,  Section  9.1.1.     Provided  that  the  relation 
between  load  and  load  effect  is  linear,  or  nearly  so,   it  makes  no  difference  when  the 
load  factors  are  applied.  Hov/ever,  for  certain  non-linear  problems,  such  as  stability 
checks,  it  is  unconservative  to  factor  the  load  effect. 

5.1  Overall  Resistance  Factors 

For  consistency  with  Galambos'   study,  criteria  taking  the  form  of  Equation  5.6  shall 

be  referred  to  as  load  and  resistance  factor  design  (LRFD) .     In  this  regard,  it  might  be 

noted  that  the  coefficients  appended  to  terms  R^and        in  Equations  2.6a  and  2,6b  are 

tantamount  to  resistance  and  load  factors  when  the  criterion  takes  the  form  of  Equation 

5.6.     Thus,  the  reliability-based  design  problem  is  to  find  an  appropriate  mean  resistance 

and  load  and  resistance  factors  for  a  prescribed  set  of    Q  ,  Q      ,6  and  load  ratios. 

R       S  . 
1 

The  same  c.o.v.   that  were  specified  for  calibration  in  the  previous  section  are  also  used 
here. 

The  resistance  and  load  factors  for  the  gravity  load  combination  (D  +  L)  are  shown 
as  functions  of  L/D  in  Figure  5.1  for  flexure,  Figure  5.2  for  shear,  and  Figure  5.3  for 
thrust  and  bending.     Since  the  calibration  showed  that  6  was  not  constant  for  existing 
designs,  the  sensitivity  of  ((>,        and        to  3  is  also  examined  by  showing  their  variation 
for  two  different  but  reasonable  values  of  3  for  the  flexure  case. 

Several  points  are  worth  noting  from  Figures  5.1  through  5.3.     The  most  important  is 

that  the  load  factors        ^rid  y     are  quite  insensitive  to  Q  ,  as  may  be  seen  by  comparing 

D  L  K 

the  flexure  and  shear  cases.     A  certain  amount  of  coupling  between  the  resistance  and 

load  factors  exists  because  the  structural  safety  condition  is  described  by  one  inequality 

out  of  which  (at  most)  one  independent  piece  of  information  (overall  safety  factor)  can 

be  extracted.     The  fact  that  this  coupling  appears  relatively  weak  will  have  some  important 

implications  in  developing  Level  I  design  criteria  to  be  discussed  subsequently.  The 

increase  in  y^  with  increasing  L/D  occurs  because  as  L/D  increases  the  live  load  with  its 

associated  higher  variability  comprises  a  larger  portion  of  the  total  design  load.  The 
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Figure  5.1  -  Load  and  Resistance  Factors 
for  Flexure  (U  =  D  +  L) 
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Figure  5.2  -  Load  and  Resistance  Factors  for 
Shear  (U  =  D  +  L) 
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Figure  5.3  -  Load  and  Resistance  Factors  for 

Compression  and  Bending  (U  =  D  +  L) 
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fact  that  Yj^  depends  on  L/D  when  3  is  prescribed  also  means  that  in  existing  design 
criteria  as  well  as  Level  I  methods,  there  will  be  some  violation  of  the  reliability 
requirement  for  certain  values  of  L/D  (unless,  of  course,        is  chosen  as  the  asymptotic 
value  at  large  L/D  ;   in  this  case,  however,  the  designs  at  other  L/D  would  be  grossly 
conservative) . 

It  might  be  observed  from  Figures  5.1  through  5.3  that        and  (fi  appear  to  be  virtually 
independent  of  L/D  .     Moreover,  the  appropriate  dead  load  factor  ranges  from  about  1.04 
to  1.12,  and  since  D'   ~  D,  the  load  factor  against  nominal  dead  load  should  be  about  1.1. 
This  is  because        is  quite  small  compared  to  other  variabilities  in  the  design  problem. 
However,  existing  dead  load  factors  range  from  1.4  to  1.7,  which  may  result  in  design 
penalties  when  L/D  is  small.     It  is  also  apparent  that  the  resistance  factors  are  con- 
siderably lower  than  those  in  existing  standards.     This  is  true  for  nominal  c|)-factors  as 
well,  as  may  be  seen  by  multiplying  the  (t-values  in  Figures  5.1  -  5.3  by  the  appropriate 
R/R'. 


This  analysis  suggests  that  current  load  factors  and  resistance  factors  tend  to  be 

too  high.     Although  the  combination  of  the  two  factors  may  give  results  that  are  about 

right,  the  current  tendency  to  lump  most  of  the  effect  of  design  uncertainty  on  the  load 
side  of  the  design  equation  clearly  is  unjustified. 

Resistance  and  load  factors  for  the  gravity  plus  wind  load  case  (D  +  L^  +  W)  are  shown 

in  Figure  5.4  for  various  L^/D  and  W/D  for  the  flexure  case;  the  discussion  which  follows 

also  applies  if  shear  or  compression  and  bending  are  considered  instead.     The  behavior  of  <^ 

and        is  virtually  the  same  as  when  only  gravity  loads  act  on  the  structure.     The  very 

rapid  increase  in  y      in  Figure  5.4a  occurs  because  Q      is  quite  large  (Q,      -  0.54  in  this 
I  ^I  ^I 

2  2 

case,  corresponding  to  beams  supporting  tributary  areas  of  greater  than  750  ft.     (69.7m  )). 
As  expected,  y^  also  increases  with  W/D;  this  is  also  shown  in  Figure  5.4b,  where  y^  is 
plotted  for  L^/D  =  0.25  and  0.50,  which  values,  according  to  the  load  ratio  analysis 
discussed  in  Section  4;  bracket  most  values  of  L^/D.     A  comparison  of  Figures  5.4b  and  5.1 
shows  that  y^  is  almost  the  same  as  y^  in  the  gravity  load  case;  the  reason  for  this  is 
that  fi^  and  fi^  are  fairly  close  in  value.     This  property  will  be  taken  advantage  of  when 
Level  I  load  and  resistance  factor  design  criteria  are  developed. 
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Figure  5.4  -  Load  and  Resistance  Factors  for  Flexure  (U  =  D  +  L  +  W) 
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Finally,  resistance  and  load  factors  for  the  case  where  wind  and  gravity  loads  counteract 

one  another,  W  -  (D  +  L) ,  are  shown  in  Figure  5.5  for  compression  and  bending.     The  behavior 

of  the  load  factors  is  quite  different  than  before,  both  y    and  y      being  less  than  unity, 

D 

while  Y^^  is  significantly  higher  (cf.  Figures  5.3  and  5.5).     Figure  5.5  shows  that  in 
contrast  with  present  practice,  a  portion  of  the  live  load  should  be  subtracted  along  with 
the  dead  load  from  y    W  when  computing  U* .     However,  the  decrease  in  y      with  increasing 

I 

L^/D  indicates  that  as  the  live  load  becomes  a  larger  portion    of  the  total  gravity  load, 

becomes  smaller  and  less  is  subtracted  from  y,,W. 
I  W 

It  may  also  be  observed  that  the  (fj-factors  are  considerably  higher  than  in  the  two 
previous  cases  examined,  and  appear  to  be  more  sensitive  to  W/D  (cf.  Figures  5.5  and 
5.3).     For  example,  here  (j)  appears  to  be  about  0.75  for  compression  as  opposed  to  the 
values  of  about  0.6  previously  observed.     This  is  a  potential  source  of  difficulty, 
assuming  that  it  is  desirable  for  (j)  values  to  be  the  same  for  all  load  combinations  in 
Level  I  design. 


It  is  apparent  from  these  three  cases  that  the  reliability  criteria  3  =        is  impossible 
to  satisfy  for  all  load  ratios  (L^/D  ,  W/D)  with  only  one  set  of  load  factors.  Thus, 
Level  I  formats  necessarily  result  in  some  local  violation  of  the  reliability  criterion. 


5.2  Partial  Resistance  Factors 

In  the  partial  factors  approach,  the  factored  ultimate  resistance  is  given  by  Equation 
5.3.     It  is  assumed  that  capacities  derived  from  ultimate  strength  principles  define  the 
resistance  function  required  by  Equations  2.5  and  2.8,  when  corrected  for  equation  bias. 
Thus,  for  example,  R*  for  flexure  would  be, 

R*  =  B*   (A  f  )*d*  (1  -  0.59q*)  (5.7) 

s  y 

and  similarly  for  shear.  Difficulties  arise  in  applying  Equations  2.5  through  2.8  to  the 
combined  thrust  and  bending  case,  however,  because  the  expression  for  resistance  does  not 
exist  in  closed  form.     The  derivatives  in  Equation  2.8  must  then  be  evaluated  numerically, 
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1.75 


and  while  this  presents  no  conceptual  problem,   in  practice  it  is  very  awkward.  For 
computing  partial  factors  for  columns,  then,  the  case  of  pure  thrust  will  be  considered, 
where  a  simple  closed  form  expression  for  ultimate  capacity  is  available.     As  the  load 
eccentricity  increases,  these  PF  should  approach  those  for  flexure,  provided  that  6 
remains  constant. 

The  inclusion  of  the  bias  term  B*,   is  important.     It  is  essential  that  the  total 
variability  associated  with  resistance  in  the  limit-state,  given  by  Equations  2.1a  and 
2.5,  be  the  same  regardless  of  whether  the  RF  or  PF  approach  is  selected  for  evaluating 
R*.     B,  and  Q    are  incorporated  in  the  R  and        used  to  derive  tf)  when  the  RF  approach  is 

D  R 

used.     The  most  straightforward  way  to  incorporate  their  effect  when  using  the  PF  method 
is  to  include  B  as  one  of  the  basic  design  variables.     While  this  may  appear  unusual  at 
first,  it  should  not  create  any  difficulties  in  practice,  as  all  that  would  be  required 
would  be  to  multiply  the  calculated  ultimate  strength,  e.g.,  M*,  by  the  numerical  constant 

Partial  factors  are  developed  in  the  following  sections  considering  that  the  factored 

ultimate  load  has  the  form  U*  =  y^^D  +  Yj^L-     The  PF  values  and  their  behavior  are  virtually 

the  same  when  wind  load  is  included.     As  before,  3  =  3.0  for  flexure,  6  =  3.0  for  shear 

and  3  =  3.5  for  thrust  in  tied  columns  and  it  is  emphasized  that  the  individual  c.o.v. 

are  the  same  as  those  used  to  evaluate  U  .     This  should  permit  a  direct  comparison  later 

R 

of  required  resistances  obtained  with  the  RF  and  PF  methods. 

Specific  partial  factors  derived  by  considering  the  bending  limit  state  are  shown  in 

Figure  5.6  as  functions  of  L/U  and  in  Figure  5.7  as  functions  of  pf^/f^  when  L/D  =  0.75. 

While  4)  ,  <b    and  6     appear  relatively  insenstive  to  L/D,  d    and  pf  /f  ,  A"  is  extremely 
B      y  e  "  ycu 

sensitive  to  d  .     This  occurs  because  fi,  decreases  significantly  as  d     increases.  Thus, 
n  d  n 

for  an  8  in  (203  mm)   thick  beam  or  slab  where  d^  =  6.3  in  (160  mm)        =  0.09  and  - 

7/8"  (Figure  5.7),  and  for  safety  checking  a  value  d*  =  d^  -  7/8"  (d^  -  22  mm)  should  be 

used.     On  the  other  hand,  if  b     -  20  in  (  508  mm)  and  d     -  16  in  (406  mm)  Q    =  0.04  and  d* 

n  n  a 

=  d    -  3/8"  (d    -  9.5  mm)  would  be  appropriate.     The  value  A"  does  not  appear  to  be 
n  n  " 

especially  sensitive  to  pf^/f^. 
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Figure  5.6  -  Variation  in  Partial  Factors  with  L/D  (pf^/f^  =  0.15;  6.3  in  <  d    <  16  in) 
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The  slight  decrease  in  (j)^  observed  in  Figure  5.7  reflects  the  increased  importance  of 
in  determining  flexural  capacity  for  heavily  reinforced  sections.     In  other  words,  as 
the  importance  of  a  particular  variable  increases,  its  partial  factor  tends  to  become 
more  conservative.     This  may  also  be  observed  with  d.     The  reason  that  <j)^  is  greater  than 
(j)^  is  that  in  a  properly  designed  underreinf orced  flexural  member,  f^  is  of  secondary 
importance  to  A^f    and  d  in  determining  ultimate  strength,  as  a  few  simple  calculations 
show.     On  the  other  hand,  small  variations  in  effective  depth  are  more  important  in  the 
safety  of  thin  members  than  deep  ones,  and  the  importance  of  good  workmanship  in  placing 
reinforcement  in  shallow  members  can  hardly  be  overemphasized. 

Partial  factors  derived  from  the  shear  limit  state  equation  are  shown  in  Figure  5.8 

as  functions  of  the  parameter  p  f  /v  .     With  the  exception  of  A",  the  PF  were  found  to  be 
^  V  y    c  ^  d 

insensitive  to  L/D  and  to  d  .     As  with  flexure.  A"  is  dependent  on  d     for  the  reasons  noted 

n  d  n 

earlier.     The  increase  in  6    and  decrease  in  A     observed  when  p  f  /v    increases  reflects  the 

c  y  V  y  c 

added  importance  of  steel  yield  stress  and  lessening  contribution  of  v    when  the  web  is 

c 

heavily  reinforced.     The  fact  that  (f)     generally  is  lower  than  for  flexure  while  ())  generally 

c  y 

is  higher  reflects  the  fact  that,  relatively  speaking,  concrete  strength  is  more  important 
in  shear  than  in  flexure. 


Finally,  partial  factors  appropriate  for  thrust  are  shown  in  Figure  5.9  for  different 

p  f  /0.85f  ,  where  p     =  (A'     +  A')bh.     Variability  in  b  and  h  are  incorporated  here  because 
g  y  c  g  s  S 

these  tend  to  be  more  important  for  thrust  than  in  flexure  and  shear.     The  PF  for  thrust 

are  insensitive  to  L/D  and  to  nominal  column  dimensions.     The  tolerance  A"  depends  to  a 

h 

certain  extent  on  h,  as  fi,    is  related  to  the  size  of  the  column.     Since  few  columns  would 

h 

be  less  than  12  in.    (305  ram)  in  size,   it  would  be  a  simplifying  assumption  to  set  h*  = 

h    -  1/4"  (h    -  6.4  mm)  for  safety  checking  columns.     6     increases  as  p  f  /0.85f  becomes 
n  n  c  g  y  c 

large  and  the  relative  amount  of  the  axial  load  carried  by  the  reinforcement  increases.  In 
all  cases,        is  much  less  than  for  either  flexure  or  shear,  being  in  the  range  0.6  -  0.7, 
because  of  the  relative  importance  of  f^  in  determining  axial  capacity. 


A  comparison  of  and  computed  using  RF  and  PF  is  shown  in  Figures  5.10a  and  b  for 
flexure  and  thrust,  respectively.     While  small  differences  naturally  exist  (due  in  part  to 
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Figure  5.8  -  Partial  Factors  for  Shear  (d^  =  16  in,  L/D  =  0.75) 
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Figure  5.9  -  Partial  Factors  for  Thrust  (b    =  h    =20  in,  L/D  =  0.75) 
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Using  Partial  and  Total  Resistance  Factors 
in  Limit  State  Equation 
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rounding  of  uncertainties  to  two  significant  figures  in  RF)  the  load  factors  computed  using 
the  two  approaches  are  virtually  the  same.     (This  is  also  true  for  y^,  as  in  Equation  5.4). 
It  may  therefore  be  concluded  that  the  choice  of  whether  to  use  RF  or  PF  in  design  codes 
may  be  deferred,  while  the  load  factors  derived  from  either  method  may  be  implemented 
without  delay.     This  is  encouraging  in  view  of  the  current  lively  discussion  on  whether  to 
use  the  RF  or  PF  approach. 

Designs  obtained  using  the  PF  and  RF  are  compared  in  Tables  5.1a,  b  and  c  for  flexure, 
shear,  and  thrust. 


Table  5.1  -  Comparison  of  RF  and  PF  Designs. 


(5.1a)  Flexure:     Design  Parameter  B  A  f  d  /D    when  pf  /f    =  0.15 

s  y  n  y  c 


L/D 

0. 

25 

0.50 

0.75 

1 

00 

2 

00 

cov(R) 

=  0.13 

RF 

2. 

199 

2.691 

3.214 

3 

753 

5 

961 

PF 

2. 

195 

2.687 

3.209 

3 

747 

5 

951 

cov(R) 

=  0.16 

RF 

2. 

376 

2.903 

3.459 

4 

031 

6 

372 

PF 

2. 

394 

2.922 

3.479 

4 

052 

6 

399 

(5.1b)  Shear:  Design  Parameter  B  v  bd  /D  when  p  f  /v  =  1.0 
  c    n  v  y  c 


L/D 

0.25 

0.50 

0.  75 

1.00 

2.00 

RF(fi„  =  0.22) 
K 

PF 

1.276 
1.177 

1.554 
1.435 

1.845 
1.  706 

2.143 
1.985 

3.362 
3.126 

(5.1c)  Thrust:     Design  Parameter  B  0.85f  bh/D  when  p  f  /0.85f    =  0.35 

c  g  y  c 


L/D 

0.25 

0.50 

0.75 

1.0 

2.0 

RF  (f^  =0.17) 

R 

PF 

1.812 
1.739 

2.216 
2.131 

2.568 
2.547 

3.088 
2.977 

4.903 
4.740 
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The  small  variation  in        for  each  limit  state  discussed  in  Section  3.2  is  sufficient  to  cause 
the  differences  between  the  RF  and  PF  methods  observed  in  Table  5.1.     However,  since  these 
differences  seldom  exceed  a  few  percent,  it  may  be  concluded  that  the  two  methods  would 
result  in  essentially  the  same  designs  under  gravity  loading  and  that  the  PF  format  would 
not  penalize  reinforced  concrete  designs. 

A  primary  advantage  to  the  PF  approach,  at  least  in  a  conceptual  sense,  is  that  the 
safety  factors  are  applied  directly  to  the  individual  design  variables  whose  randomness  is 
the  source  of  variability  in  resistance.     Presumably,  if  one  variable  was  particularly 
important  to  safety  its  PF  would  reflect  that  importance. 

It  has  been  proposed  [12,26]   that  this  method  be  implemented  in  design  by  prescribing 
one  set  of  factors,  (j)^,(j)^,A^',A|j,  etc.,  that  would  be  used  for  all  ultimate  limit  states,  and 
that  this  would  simplify  design.     Although  the  importance  of  including  equation  bias  and 
variability  arising  therefrom  has  not  been  specifically  mentioned,  its  inclusion  presents 
no  particular  problem  as  long  as  B  is  accepted  as  one  of  the  basic  design  variables.  As 
was  the  case  with  RF  design,  however,  PF  consistent  with  a  specified  reliability  differ 
for  the  various  limit  states.     This  is  due  to  differences  in  the  relative  importance  of 
steel  and  concrete  strength  and  member  geometry  in  each  limit  state,  as  well  as  the 
different  6  required  in  each  case.     The  use  of  the  same  6  for  all  limit  states  reduces, 
but  does  not  entirely  eliminate,  this  difference.     Thus,  the  use  of  one  set  of  PF  for 
flexure,  shear  and  thrust  would  result  in  local  violations  of  the  reliability 

constraint.     For  example,  if  (J)^  =  0.85  -  0.9,  as  suggested  by  the  flexure  and  shear  solutions, 

the  corresponding  safety  index  for  thrust  would  be  about  3.0  rather  than  the  3.5  required. 

Similarly,  the  use  of  one  bar  placement  tolerance,  say  A"  =  3/8"  (9.5  mm)  would  be  unconser- 

d 

vative  for  shallow  beams  and  their  slabs.  On  the  other  hand,  the  effect  of  small  changes  in 
individual  variables  on  the  total  variability  in  resistance  tends  to  get  damped  out  in  the 
RF  approach.     This  is  particularly  true  for  shear,  where  variability  associated  with  the 
prediction  equation  itself  is  a  significant  component  of  the  total  variability. 


Although  the  partial  factors  procedure  will  not  be  considered  further  for  the  above 
reasons,  with  the  exception  of  shallow  members   (h  <  8",   say)   the  following  PF  on  mean 
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values  appear  appropriate:     i>_^  =  0.84,  (j)^  =  0.65,        =  1/2"  (13  mm),        =  1/4"  (6  mm);  A^ 
=  1/8"  (3  mm).     The  factored  resistance  R*  should  include  the  term  B*  as  noted  previously. 
The  load  factors  should  be  the  same  as  for  LRFD. 


5.3  Single  Load  Factors 

As  a  final  design  format  example,  we  consider  the  equation, 

<})R    >    Y  (D  +  L)  (5.8) 

which  is  the  simplest  format  in  which  it  is  possible  to  treat  the  load  and  resistance  sides 

of  the  design  equation  separately.     For  illustration,  flexural  resistance  will  be  considered, 

wherein  Q    =0.16  and  6  =  3.0;  similar  results  are  obtained  with  shear  or  thrust.  The 
R 

c.o.v.  in  the  total  load  is. 


Vl  =    I^b  a/D)^.J]^/2  (5.9) 


and        being  the  same  as  before. 


The  behavior  of      and  y  is  shown  in  Figure  5.11,  along  with  the        and  y^  detemined 
as   before.     As  would  be  expected,  y  falls  between  y„  and  y  ,  approaching  y  at  small  L/D. 

JJ  Li 

The  resistance  factor  (j)  is  exactly  the  same  as  the  LRFD  value.     This  indicates  again  that 
the  functional  form  of  the  resistance  or  load  side  of  the  design  equation  may  be  altered 
without  significantly  affecting  the  safety  factors  on  variables  on  the  opposite  side.  It 
may  also  be  recalled  from  the  previous  section  that  y^  and  y^  were  the  same  for  both  PF  and 
overall  <j)-factor  approaches.     This  suggests  that  standard-writing  groups  concerned  with 
resistances  and  those  concerned  with  loads  could  proceed  independently  in  developing  criteria 
without  agreeing  precisely  on  an  overall  format  for  design  (e.g.  PF  vs.  RF) .     It  is  only 
necessary  that  they  agree  on  the  underlying  measures  of  uncertainty        and  the  reliability 
analysis  methods  through  which        and  6  are  incorporated  to  derive  design  criteria. 
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It  was  observed  in  Table  4.1,  that  when  only  one  (constant)  load  factor  is  applied  to 
the  sian  of  load  effects  in  a  Level  I  procedure,  as  in  Equation  5.4,  the  associated  value  of 
3  varies  more  with  L/D  than  if  the  load  factors  are  separated,  as  in  Equation  5.5.  The 
Level  I  procedure,   to  be  discussed  subsequently,  requires  that  a  set  of  load  and  resistance 
factors  be  chosen  that  minimize  the  reliability  constraint  violation,  and  this  may  be 
better  achieved  with  the  separate  load  factor  procedure.     Additional  practical  difficulties 
with  Equation  5.4  arise  when  lateral  forces  counteract  gravity  loads,  as  the  c.o.v.  in  (W  - 
D  -  L)  may  become  very  large  in  this  case.     If  one  overall  y  is  desired,  it  appears,  on  the 
basis  of  Figure  5.11  and  supplementary  calculations  considering  D  +  L  +  W,  that  y  should  be 
in  the  range  1.25  -  1.3.     It  seems  worthwhile  to  re-emphasize  the  conceptual  advantages, 
discussed  at  the  beginning  of  Section  5,  that  this  approach  has  for  Level  II  design  procedur 
wherein  variation  of  <p  and  y  with  (L/D,  W/D,   '■')  would  be  an  acceptable  feature. 
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6.  PRACTICAL  LOAD  AND  RESISTANCE  FACTOR  DESIGN 


At  present,  Level  II  methods  are  useful  primarily  to  technical  committees  charged  with 
developing  design  standards.     However,   for  simplicity  and  continuity  with  existing  formats, 
it  is  necessary  to  eliminate  the  dependence  of  the  safety  factors  on  the  load  ratios  and  to 
prescribe  a  minimal  number  of  such  factors  suitable  for  use  in  everyday  design.     In  the 
short  term,   the  LRFD  approach  appears  most  appropriate  for  Level  I  formats  in  the  United 
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States  because  of  its  consistency  with  the  existing  ACI  [10]  and  the  proposed  AISC  provisions 
[60]. 


The  maximum  total  load  on  a  structure  or  component  during  its  life  may  result  from  any 
one  of  many  different  possible  load  combinations.     It  has  been  found   [63]  that  since  the 
probability  of  a  joint  occurrence  of  two  or  more  extreme  loads  is  negligible,  the  maximum 
total  load  effect  results  when  one  variable  load  assumes  its  maximum  lifetime  value  while 
others  are  at  their  instantaneous  values.     The  factored  ultimate  load  U*  thus  becomes, 

n  _ 

U*  =        D  +  Yq  Q  +  S.  (6.1) 

i  i 

where  D  is  mean  permanent  (dead)  load,  Q  is  the  mean  maximum  variable  load  during  some 
suitable  period,  e.g.,  the  life  of  the  structure,  and        is  the  mean  instantaneous  or 
short-term  load  which  acts  simultaneously  with  Q.     In  the  general  case,  the  individual 
loads  under  consideration  must  be  rotated  in  the  load  equation,  each  load  taking  the  position 
of  the  maximum  mean  variable  load  Q,  while  the  remaining  loads  are  assigned  the  value  S^. 
So  that  a  direct  comparison  can  be  made  with  the  requirements  of  Sections  9.2  and  9.3  of 
ACI  318-71,   it  will  be  assumed  that  Q  equals  either  L  or  W,  and  that  S^  =  L^.  This  by  no 
means  encompasses  all  possible  load  combinations,  but  will  enable  the  concepts  required  for 
Level  I  criteria  development  to  be  demonstrated.     The  factored  ultimate  resistance  R* 
against  which  U*  is  compared  is  simply  R*  =  (})R,  in  accordance  with  Equation  5.2. 


The  variation  in  (<i) ,  y  ,  y  ,  y      )  with  the  load  ratios  in  Section  5,  shows  that  if  a 

D  Q 

constant  set  of  these  values  are  prescribed,  the  reliability  constraint  will  be  violated 

locally.     While  such  violations  are  unavoidable  in  Level  I  design,  it  nevertheless  should 

be  possible  to  select  one  set  of  LRF  that  minimizes  the  extent  of  this  violation  when 

considered  over  all  likely  load  ratio  combinations.     One  way  to  do  this  is  to  select  a  set 

2 

of  LRF  that  minimize  the  quantity  EEp..(6..-3  )   ,  where  6..   is  the  safety  index  obtained 

using  a  prescribed  LRF  set  for  load  combination  ij  occurring  with  probability  p     ,  and  6 

ij  o 

is  the  target  reliability  [58].     However,  this  approach  is  numerically  difficult  to  apply 
with  the  reliability  procedure  given  by  Eqs.  2.5-2.8  because  their  solutions  are  not  available 
in  closed  form. 
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As  an  alternative  procedure,  a  function  of  the  form 


I    -  [R     (3  )  -  R  (Safety  factors;  load  ratios  )]' 

■■-J         -LJ         O  j[j 


(6.2) 


is  minimized.     In  Eq.   6.2,  R. ^ (6^)  is  the  mean  resistance  required  according  to  Eqs.  2.5- 
2.8  for  a  target  reliability  6^  and  the  load  combination  ij .     The  function  R  is  the  mean 
resistance  furnished  by  the  Level  I  design  procedure,  e.g.,   if  LRFD  is  used,  R  =  (y^^D  + 
Y  L  +  Y  W 

L         W  Note  that  this  minimization  procedure,  as  well  as  the  one  described  in  the 

preceding  paragraph,  requires  the  selection,  a  priori,  of  a  Level  I  design  format.  Indeed, 

the  specific  safety  factors  extracted  will  depend  on  the  functional  form  of  R  and  constraints 

placed  upon  it.     Compared  to  the  procedure  outlined  in  the  previous  paragraph,  however,  the 

minimization  of  Eq.   6.2  is  relatively  simple,  since  R. . (6  )  are  tabulated  during  the  Level 

II  development  discussed  in  Section  5.     Moreover,  it  was  found  by  computing  ZZp     (g  - 

2  ij  ij 

e^)     for  small  perturbations  in  the  LRF  values  obtained  from  Equation  6.2  that  the  LRF  that 

minimize  both  expressions  are  within  a  few  percent  of  one  another. 

The  determination  of  an  appropriate  set  of  mean  load  ratios  for  this  minimization  and 
the  assignment  of  relative  likelihoods  p_  to  them  cannot  be  done  objectively  given  the 
present  state  of  knowledge.     In  a  recent  study  [58],  a  set  for  all  limit  states  and  materials 
was  selected  subjectively.     However,  since  the  present  study  is  concerned  with  reinforced 
concrete  structures,   the  L/D  and  W/D  ratios  would  be  smaller  than  those  previously  reported. 
Accordingly,  the  values  listed  in  Table  6.1a  were  assigned  for  the  (D  +  L)  case; 


Table  6.1a  -  Mean  Load  Ratios  and  Weights  for  D  +  L 


L/D 

0.25 

0.  50 

0.75 

1.0 

2.0 

Pi 

0.22 

0.31 

0.25 

0.14 

0.08 

and  when  live  load  is  combined  with  wind,  as  in  Table  6.1b; 


Table  6.1b  -  Mean  Load  Ratios  and  Weights  for  D  +  L^  +  W 


L^/D 

0.10 

0.25 

0.50 

0.75 

1.00 

Pi 

0.21 

0.44 

0.29 

0.04 

0.02 
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For  each  L^/D,  these  p.  must  be  apportioned  to  possible  W/D  values  also.     On  the  basis  of  an 
analysis  of  a  multistory  moment  resisting  reinforced  concrete  frame,  using  loads  defined  in 
Appendix  A,  equal  weights  p.,  were  assigned  all  W/D  values  between  0.5  and  2.5.     A  study 
was  also  conducted  to  determine  the  sensitivity  of  the  Level  I  criteria  to  the  load  ratio 
and  probability  assignment.     It  was  found  that  although  the  load  and  resistance  factors 
depend  on  the  range  of  L/D,  W/D  assumed,  surprisingly  enough  they  do  not  appear  to  be 
especially  sensitive  to  the  assignment  of  p^_.  within  that  range. 

Using  the  LRFD  format  for  load  combination  k,  Eq.  6.2  has  the  form  (some  factors  may 
be  zero) , 

Considering  more  than  one  load  combination  simultaneously,  the  function  to  be  minimized  is. 


in  which  (<j) ,Yp,>Y^>Y-r    )  are  the  required  Level  I  load  and  resistance  factor  set. 
D    Q    L  J 

Certain  additional  restrictions  are  placed  on  the  minimization  problem.     First  of  all, 
it  was  apparent  from  the  results  in  Section  5  that        =  1.10  in  all  cases  where  gravity  and 
lateral  loads  are  additive.     Secondly,  it  was  observed  (cf.  Figures  5.1  and  5.4b)  that  the 
load  factors  applied  to  the  maximum  live  load  and  wind  loads  were  quite  similar  because 
n    and  Q.    are  close  in  magnitude.     Thus,  it  appears  reasonable  to  set  y     (=  Yt   =  Y,,)  constant 
in  Equation  6.1.     In  other  words,  the  load  factor  applied  to  the  maximvm  variable  load  in 
the  load  equation  should  be  the  same,  regardless  of  the  load  combination.  Intuitively, 
this  appears  reasonable  because  the  c.o.v.  associated  with  most  maximum  variable  loads 
(with  the  exception  of  earthquake  loads)  is  generally  about  0.30.     However,   this  represents 
a  slight  departure  from  current  ACI  practice,  wherein  according  to  Eq.  9-1  and  9-2,  Yt  - 
1.7  and  y^^  =  1.  275,  respectively. 

Minimization  of  Equations  6.3  is  first  performed  separately  for  flexure,  shear,  and 
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axial  compression  and  bending,  considering  the  two  load  combinations  (EH-L,  D+W+L^)  wViere  load 
effects  are  additive.     The  results  of  this  analysis  are  presented  in  the  first  four  columns 
of  Table  6.2. 


Table  6.2  -     Optimal  Load  and  Resistance  Factors 


Limit  State 

(Yq  =  1.55;  y^^  =  1.6) 

0.77 

1.10 

1 

59 

1.63 

0.76 

0.71 

1.10 

1 

54 

1.59 

0.71 

Shear 

0.59 

1.10 

1 

47 

1.50 

0.61 

Compression  &  Bending 
Tied  Columns 

0.63 

1.10 

1 

59 

1.63 

0.62 

Spiral  Columns 

0.67 

1.10 

1 

55 

1.59 

0.66 

Selecting  y  =  1.55  and  y  =1.6,  thus,  the  values  of  <})  which  minimize  Equations  6.3  are 
listed  in  the  last  column  of  Table  6.1. 


In  order  to  determine  an  appropriate  set  of  load  factors  for  the  third  load  combination, 

when  gravity  and  wind  forces  counteract  one  another  (W-D-L^),  we  require  only  that  the  (i>- 

values  be  identical  to  those  listed  in  Table  6.1.     No  other  constraints  are  the  values  of 

Y  ,Yy.    and  Yj       are  placed  (the  "minus"  superscript  means  that  the  associated  loads  are  to 
X    D  L  J 

be  subtracted  in  the  load  combination).     Minimizing  Equation  6.3  loads  to  Yq  ~  1.60,  y^ 

=  0.8  and  y^^  =  0.60.     It  would  be  conservative  to  neglect  the  latter  term.     The  load  criteria 

analogous  to  Equations  9-1,  9-2  and  9-3  in  ACI  318-71  would  thus  be, 

U*  1.1  D  +  1.55  L  (6.4a) 
U*  =  1.1  D  +  1.6  Lj  +  1.55W  (6.4b) 
U*  =  1.60  W  -  (0.8  D  +  0.6  L^)  .(6.4c) 


The  <{)  values  in  the  last  column  of  Table  6.1  and  y  in  Equations  6.4  are  associated  with 
mean  resistance  and  loads.     To  obtain  appropriate  values  to  be  applied  to  nominal  quantities. 


simply  multiply  <i>  by  R/R'  and  y^  by  ^-j^/^^- 
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The  typical  violation  of  the  reliability  constraint  arising  from  the  selecting  of  a 
Level  I  criteria  is  illustrated  in  Table  6.3,  where  3,.  for  the  shear  criterion  0.61R  = 
I.ID  +  1.55W  +  1.6L^  are  given.     The  target  value  3^  was  3.0. 

Table  6.3    -  3  Achieved  by  Level  I  Criterion  for  Shear 


L^/D 

W/D 

0.10 

0.25 

0.50 

0 

75 

1.0 

0.5 

2.98 

3.02 

2.98 

2 

90 

2.82 

1.0 

3.01 

3.06 

3.07 

3 

03 

2.98 

1.5 

3.00 

3.05 

3.09 

3 

08 

3.05 

2.5 

2.97 

3.01 

3.06 

3 

08 

3.09 

Using  the  p..  in  Table  6.1,  the  mean  and  variance  of  the  safety  index  achieved  are 
listed  in  Table  6.4. 


Table  6.4  -  Variability  in  3  Caused  by  Level  I  Criteria 


Limit  State 

3 

o 

Load  Combination 

E[3]=e 

E[(3-0)^] 

Flexure 
nj^=0.16(0.13) 

3.0 

D  +  L 

D  +  L  +  W 

2.99(3.03) 
3 .01(3,01) 

0,00975(0,00974) 
0,00322(0,00536) 

Shear 

3.0 

D  +  L 

D  +  L  +  W 

2  ,96 

3  ,03 

0,00876 
0,00174 

Compression 
and  Bending 
(Tied  Columns) 

Compression 
and  Bending 
(Spiral 
Columns) 

3.5 
3.2 

D  +  L 

D  +  L  +  W 

D  +  L 

D  +  L  +  W 

3.55 
3  ,54 
3  ,23 
3  .24 

0,00698 
0,00408 
0,00838 
0.00327 

It  may  be  observed  that  3  is  within  0.1  of  3^  in  all  cases.    This  may  be  contrasted  with 
the  scatter  in  3  values  associated  with  existing  design  criteria,  discussed  in  Chapter  4. 

It  is  worth  reemphasizing  that,  in  accordance  with  Equation  6.2,  the  Level  I  load  and 
resistance  factors  depend  on  the  format  that  is  assumed,  a  priori,  for  the  design  equation. 
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In  principle,  this  could  range  from  a  single  equation  with  one  overall  safety  factor  to  the 

complex  formats  currently  under  consideration  by  CEB  [12]     and  other  European  code-writing 

groups.     As  with  multiple  regression  analysis,  the  more  independent  factors  that  are 

assigned,  the  less  the  reliability  constraint  will  be  violated,   i.e.,  b"  in  Table  6.4  will 
—  2 

approach  3^  and  E[(3-3)   ]  will  approach  zero.     This  militates  against  the  use  of  one 
overall  safety  factor  or  load  factor  in  Level  I  criteria  development,  and  would  suggest 
that  as  many  factors  be  assigned  as  would  be  practical  and  easy  to  handle  for  everyday 
design  usage. 
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7.  CONCLUSIONS 


Reliability  analysis  methods  have  been  employed  in  this  study  to  facilitate  the  selection 
of  criteria  for  reinforced  concrete  design.     These  methods,  which  are  based  on  probability 
theory,   provide  a  logical  basis  for  determining  the  manner  in  which  uncertainties  in  resistance 
and  loads  affect  design  safety  and  how  their  effects  should  be  controlled  in  building 
codes.     Using  data  obtained  from  a  comprehensive  analysis  of  uncertainty  measures  (described 
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in  Appendix  A) ,  safety  indices  associated  with  existing  reinforced  concrete  designs  were 
evaluated.     This  provided  target  reliabilities  for  subsequent  criteria  development.  Finally, 
Level  II  design  criteria  commensurate  with  levels  of  uncertainty  and  required  reliability 
were  presented.     Simplification  of  these  led  to  practical  reliability  based  criteria  which 
retain  the  relatively  simple  characteristics  of  existing  criteria  and  yet  have  a  well 
established  and  documented  rationale. 

This  study  has  presented  a  rather  broad  treatment  of  the  use  of  reliability  methods  in 
developing  criteria  for  structural  design  against  various  ultimate  limit  states.     By  examining 
the  sensitivity  of  the  solutions  to  certain  assumptious  or  simplifications  a  decision  can 
be  made  as  to  whether  further  examination  is  worthwhile.     Additional  studies  on  the  treatment 
of  load  combinations  appears  desirable  in  the  short  term,  for  example.  Additional  data  are 
also  needed  to  determine  appropriate  load  ratios  and  relative  frequencies  as  in  Table  6.1. 
However,  since  standards  development  is  evolutionary  in  nature,  many  of  the  concepts  and 
results  discussed  herein  can  be  implemented  concurrent  to  supplementary  research. 

The  following  specific  conclusions  can  be  drawn: 

1.  Analysis  of  uncertainties  in  resistance  for  reinforced  concrete  members  shows 
that,  with  the  exception  of  shear,  they  are  about  the  same  as  for  steel  members  [60]. 
c.o.v.   in  resistance  may  approach,  and  occasionally  exceed,   those  for  loads. 

2.  Levels  of  safety  for  reinforced  concrete  designs,  measured  by  3,  average  from 
2.5  to  4,  depending  on  the  limit  state  and  tributary  loaded  area.     This  is  the  same  order 
as  found  by  Galambos  [28]   for  steel  designs.     It  is  less  than  found  by  Siu,  et.al  [58] 
because  the  ANSI  live  load  is  less  than  the  Canadian  [43]  live  load. 

3.  Reliabilities  associated  with  reinforced  concrete  beams  in  flexure  and  columns 
appear  appropriate  in  a  relative  sense.     Current  shear  provisions  appear  inadequate;  the 
nominal  <))'  =  0.85  is  insufficient  to  account  for  the  large  (0.22)  c.o.v.  in  shear  resistance. 


4.     Reliabilities  associated  with  the  current  loading  provisions 

Chapter  9  of  ACI  318-71  generally  appear  risk-consistent,  although  there 
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contained  in 

is  evidence  that 


the  reliability  against  wind  may  be  slightly  less  than  against  gravity  loads. 

5.  Load  and  resistance  factors  for  a  limit  state  which  are  consistent  with  a 
prescribed  reliability  depend  on  mean  load  ratios.     The  y     factors  do  not  appear  to  be 
strongly  dependent  on      ,  however,  suggesting  that  the  coupling  between  resistance  and  load 
sides  of  the  design  equation  is  weak. 

6.  Safety  factors  in  existing  standards  tend  to  be  allocated  improperly  to  resistance 
and  load  sides  of  the  design  equation  because  of  failure  to  account  for  relative  magnitudes 

of  uncertainty  in  resistance  and  load.     Calibration  assures  that  existing  and  reliability- 
based  designs  will  be  about  equal;  although  the  resistance  and  load  factors  in  the  latter 
will  both  tend  to  be  lower,  their  net  effect  will  be  about  the  same. 

7.  The  load  and  overall  resistance  factor  (LRFD)  format 

(4)R  >  I  (y.S') 

appears  more  attractive  for  American  designers  than  either  the  partial  resistance  factor  or 
single  load  factor  approaches.     However,  since  the  factors  obtained  on  one  side  of  the 
design  equation  are  not  sensitive  to  the  format  chosen  for  the  other  side,  provided  that 
measures  of  uncertainty  and  6  are  agreed  upon,  the  question  of  whether  to  use  overall  or 
partial  resistance  factors  need  not  be  agreed  upon  prior  to  load  criteria  development. 

8.  Level  I  criteria  depend  on  the  format  selected  and  constraints  imposed.  Practical 
criteria  can  be  developed  from  reliability  bases  which  are  no  more  complex  than  design 
procedures  currently  in  use.     As  an  example,  for  reinforced  concrete  design,  the  following 
load  criteria  and  resistance  factors  could  replace  those  in  ACI  318-71.     Note  that  the  load 
and  resistance  factors  are  applied  to  mean  values. 
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Load : 

U*=1.1D+1.55L 

U*  =  1.1D+  1.55  W+  1,6L^ 

U*  =  1.6    W  -  (0.8  D  +  0.6  L^) 

Resistance 

Flexure  <^  =  0.73 

Shear  ^  =  0.61 

Compression  and  Bending  (Tied  Columns)        (}>  =  0.62 

Compression  and  Bending  (Spiral  Columns)     <\>  =  0.66 
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APPENDIX  -  ANALYSIS  OF  UNCERTAINTIES 

The  analysis  of  uncertainties,  which  was  summarized  in  Tables  3,1  and  3.2  in  Section 

3  is  described  in  detail  in  this  Appendix.     The  estimation  of  measures  of  uncertainty  in 

such  basic  resistance  variables  as  f  ,  f  ,  etc.   is  discussed  in  section  A.l.     These  are 

y  c 

combined  to  determine  total  variability  Q  and  mean- to-nominal  ratio  R/R'   in  A. 2  for  flexure, 
shear,  and  thrust  combined  with  bending.     A  similar  evaluation  of  load  statistics  is 
presented  in  A. 3.     ACI  318-71  standard  notation  is  used  whenever  possible. 

A.l    Estimation  of  Individual  Measures  of  Uncertainty 

Dispersion  in  the  basic  variables  such  as  concrete  strength,  reinforcement  yield 
stress,  placement  of  bars,  and  member  dimensions  all  contribute  to  the  variability  in 
resistance  of  reinforced  concrete  members. 

Dispersion  in  concrete  strength  is  caused  by  variability  in  material  properties,  mix 

quality  control,  and  in  transporting,  placing  and  field  curing  the  concrete.     The  basic 

variability  6^     is  determined  from  standard  cylinder  test  data  which  relate  quality  control 
c 

to    6^    by  [50] 
c 

Good  quality 
Average  quality 
Poor  quality 

Additional  uncertainty  in  f    arises  from  unavoidable  variations  in  quality  control  at  the 

c 

batch  plant  in  preparing  the  concrete.     The  dispersion  of  the  observed  sample  mean  strengths 

f    of  concrete  delivered  during  12  successive  monthly  periods  is  about  0.07   [50].  The 
c 

design  compressive  strength  f^  is  intended  to  be  a  conservative  value,  so  that  the  probability 

of  actually  achieving  a  lesser  strength  is  small.     The  mean  cylinder  strength  for  a  batch 

of  concrete  f  ,   therefore  must  be  somewhat  higher  than  f  in  order  to  satisfy  section 

c,  cyl  c 

4.2.2  of  ACI  318.71.     For  a  nominal  4000  psi  mix  with  average  quality  control  (6^    =  0.15 

c 

or  =  600  psi)   f  ,   =  1.2  f   to  satisfy  this  requirement, 

f  ^  c,  cyl  c 


0.10  <        <  0.15 

c 

0.15  <  <5j  <  0.20 

c 

0.20  <  6, 
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It  is  widely  recognized  that  the  strength  of  concrete  in  situ,  which  would  be  the 
relevant  value  to  use  for  reliability  calculations,  is  different  from  the  standard  cylinder 
test  value.     Data  presented  by  Bloem  [8]   show  the  relation  between  in  situ  strength  (in 
slabs),  as  determined  from  drilled  cores,  and  the  standard  cylinder  test  strength.  Data 
given  in  Table  2  of  RefjSare  plotted  in  Figure  A.l.     Regression  analyses  of  the  "actual" 
drilled  core  strength  f  on  standard  cylinder  strength  f  yielded, 

O  y    C  O  IT  6  C 


{ 


0.67  f  ,  -  Poor  cure 

f                                    /  "^y^ 

c,  core             =         ^  (A. la) 

0.84  f  T  -  Good  cure 

c,  cyl 


with  a  (see  Equation  3.6)  =  0.098  and  0.071  for  poor  and  good  field  curing,  respectively. 
According  to  Equation  3.6,  then,  the  contribution  of  field  curing  to  concrete  strength 
variability  ranges  from  0.08  to  0.15,  with  an  average  value  of  about  0.12.     An  alternate 
regression  analysis  with  two  undetermined  constants  yielded 


{662  +  0. 
596  +  0. 


.49  f        T   -  Poor  cure 
f  1  c,cyl 

c,core  =  \  (A. lb) 

67  f  ,  -  Good  cure 

c,  cyl 


Within  the  range  3000  <  f^  <  4000  psi,   the  results  are  practically  the  same. 


Other  factors  influence  the  variability  in  f  .     For  example,  moisture  migration 

c 

upward  and  settlement  of  aggregate  cause  vertical  or  deep  members  to  be  stronger  at  the 
bottom  than  at  the  top.     At  the  bottom  and  top  of  columns,  the  concrete  strength  was  8 
percent  higher  and  4  percent  lower,  respectively,  than  the  average  column  strength  [19]. 
Using  either  Equations  3.2a  or  3.2b,  the  contribution  of  this  segregation  to  total  variability 
in  resistance  is  0.03.     Other  factors  are  more  difficult  to  assess;  the  sustained  load  strength 
of  concrete  is  about  85  percent  of  its  short  term  strength  for  example,  and  low  strain 
rates  may  only  yield  about  85  percent  of  the  strength  obtained  using  ASTM  test  methods 
[34,50,53].     On  the  other  hand,  concrete  gains  strength  beyond  its  28-day  design  strength; 
after  one  year,  Bloem  [8]   found  this  increase  could  be  in  excess  of  30  percent.  Moreover, 
confining  the  concrete  elevates  its  strength  also.     These  factors  tend  to  offset  one 
another . 
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Figure  A.l  -  Variability  in  Concrete  Strength 
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The  strength  of  concrete  in  situ        is  thus, 
f  ,  f 

f    =  f     ^'  (other  factors)  (A. 2a) 

c        c      r  r  . 

c  c,  cyl 

Using  Eq.  A.l  and  recalling  that  f^  1.2  f^,  the  mean  strength  f^  is 


f     for  poor  field  cure 
f^^  (  '  (A.2b) 

-  01  f^  for  good  field  cure 


The  total  in  situ  variability  is. 


=     [6^    +  (0.07)^  +  (0.12)^  +  (0.03)^]^''^  (A. 2c) 

c  c 

This  leads  to  values  of          which  range  from  0.174  to  0.207  as  6^     increases  from  0.10  to 

c  c 

0.15. 


The  variability  in  concrete  strength  in  existing  buildings  was  determined  by  Drysdale  [19] 

using  ultrasonic  techniques.     Variability  in  pulse  velocity  appears  to  correlate  quite 

well  to  that  in  f^.     For  1145  columns  in  14  buildings  Drysdale  found  that  varied  from 

c 

0.07  to  0.30.     In  four  buildings,  Q.^    was  in  the  range  0.05  to  0.10;  in  eight  0.10  to 

c 

0.20;  and  in  two,  exceeded  0.20.     Accompanying  standard  cylinder  test  data  were  unavailable. 

c 

In  an  earlier  study  using  similar  techniques,  Tso  and  Zelman  [62]   found  that  test  cylinder 
strengths  did  not  correlate  well  with  in  situ  strength,  undoubtedly  due  in  part  to  the 
additional  sources  of  variability  discussed  in  the  preceding  paragraphs. 


Tensile  strength  in  concrete  is  commonly  measured  by  the  modulus  of  rupture  or  the 

cylinder  splitting  strength.     In  one  set  of  tests   [59],  the  c.o.v,  in  the  modulus  of  rupture  ranged  froi 

0.13  -  0.18;   in  another   [47]  the  c.o.v.  in  the  splitting  strength  varied  from  0,06  to  0.14. 

These  appear  to  be  roughly  equivalent  to  the  6^     found  from  standard  cylinder  compression 

c 

tests;  additional  sources  of  uncertainty  tend  to  increase  the  variability  in  tensile 
strength  as  in  compression.  Thus,  variabilities  in  tensile  strength  and  compressive 
strength  do  not  appear  to  be  distinguishable  and  will  be  assumed  to  be  equal  in  this 
study. 
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Variability  in  reinforcement  yield  strength  is  determined  primarily  from  mill  test 

data.     Data  reported  by  Julian  [35]   for  Grade  40  reinforcement,  No.   3  through  No.  10  bars, 

—  2 
showed  that  f     =  47.7  ksi  [329  MN/m  ]   and  6^     =  0.12.     Recent  data  reported  by  Mirza  and 

MacGregor  [41]  showed  that  for  Grade  40  bars  (No,   3  through  No.  11)  f^  =  48.8  ksi  (336 

MN/m^)  and  6^     =  0.107  while  for  Grade  60  bars  (Nos.  5  through  11)  ?    =71.1  ksi  (490 

2  ^  ^ 

MN/m  )  and  6^     =  0.093.     These  data  include  an  uncertainty  due  to  size  effects.  Earlier 

y  -  2 

data  reported  by  Baker  [6  ]  showed  that  f    decreased  from  50.4  ksi  (348  MN/m  )   for  No.  3 

2 

bars  to  44.1  ksi  (304  MN/m  )   for  No.   8  bars;  6^     ranged  from  0.07  to  0.11  but  this  variation 

y 

was  not  related  to  bar  size.     On  the  basis  of  these  data,  it  appears  that  6^     ^  0.09  and 

_  y 

the  variability  in  the  mean  f^  due  to  the  size  effect  is  about  0.04  [21]. 


It  should  be  noted  that  this  6^    may  also  actually  include  the  effect  of  variability 

y 

in  the  bar  area  A  ,  since  it  is  not  clear  in  the  data  sources  cited  above  whether  f  was 
s  y 

determined  on  the  basis  of  nominal  area  A    or  actual  bar  area.     If  T  =  A  f  ,  and 

n  s  y 

A  f 

f  =  (A.  3a) 

y,  n  A^ 

2 

then  the  observed  6^  is  actually 

y,  n 

2  2  2  2 

6  S      +  &^     =  S  (A. 3b) 

L  A-  t  L 

y,  n        S  y  y 

Since  6^    is  about  0.02   [21]   this  source  of  error  should  be  negligible. 

s 

Mill  tests    upon  which  the  above  assessments  of  variability  are  made     indicate  a  bar 

yield  strength  which  is  higher  than  would  probably  be  observed  in  practice  because  of  the 

elevated  strain  rates  at  which  tests  are  made.     Julian  [35]  has  indicated  that  this  elevation 

may  be  on  the  order  of  5  -  10  percent.     This  was  borne  out  by  Allen  [2],  who  found  that 

2 

the  elevation  in  mean  value  was  about  3  ksi  (20.7  MN/m  )  for  Grade  40  reinforcement,  with 
a  c.o.v.  of  about  0.13.     If  the  yield  strength  in  the  structure  is  given  as 


f    =        f         . T  T  -  x„  (A. 4a) 

y        1    y,  mill  2 
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where        and  yl^  are  size  effect  and  strength  elevation,  respectively,  then 

2  2         _    _    _  2  2 

A      +6^    +  (x„/x,   f         .„)  A 

„2  ^1        %  ^     ^    y>  '""^^        ^2  (A.4b) 

S  (1  -  ^2/^1  \,  mill^^ 

Substituting  A      =  0.04,  6      =  0.09  and  A      =  0.13,  and  assuming  that  x    is  10  percent  of 

_  1  y  2 

X,   f  leads  to  fi^    =  0.11. 

1    y,  mill  f 

y 

Recent  data  [25]  on  variability  in  reinforced  concrete  members  dimensions  and 

bar  placement  has  shown  that  these  variabilities  have  a  definite  tendency  to  decrease  as 

the  member  size  increases.     This  is  illustrated  in  Figure  A. 2  for  member  thickness,  h.  The 

6,    reported  by  Johnson  [33]  and  Johansson  and  Warris   [32]  was  obtained  from  measurements 
h  \ 

taken  on  in  situ  slabs,  all  of  which  were  less  than  8  in  (203  mm)  thick.     The  data  reported 

by  Tso  and  Zelman  [62]  for  columns  ranging  from  12  to  30  inches  (304  to  762  mm)  shows  this 

decrease  more  clearly.     It  appears  from  Figure  A. 2  that  6,    -  0.4/h  ,  where  h    is  the 

h  n  n 

nominal  minimum  member  dimension.     Additional  confirmation  for  this  is  provided  by  Birkeland 

and  Westhoff  [7],  who  estimate  (on  the  basis  of  extensive  field  experience)  that  member 

dimensions  vary  by  4-  7/8  inch  (22  mm).     Using  Equation  3.2b  as  previously  discussed  leads 

to  6,       0.44/h  .     Finally,  Mirza  and  MacGregor   [41]  suggest  that  6,    -  0.46/h    for  slabs, 
h  n  h  n 

0.214/h    for  beams,  and  0.262/h    for  columns.     The  latter  study  also  indicates  that  the 
n  n 

variability  in  beam  width  is  quite  small. 


Additional  uncertainty  arises  from  variation  in  mean  dimensions  from  those  nominally 

specified.     The  slab  data  suggest  that  this  variability  is  about  0.04  [32,33].     Data  in 

Ref.  62,  Table  2  for  thicker  members  reveals  that  the  mean  tolerance  varies  between  -3  and 

+5  percent;  using  Equation  3.2b,   this  implies  the  variability  due  to  uncertainty  in  mean 

value  prediction  is  about  0.02;   thus,  A^  also  depends  on  the  size  of  the  member,  and 

appears  to  be  roughly  0.25/h  .     Combining  this  with  the  previous  estimate  for  6,  ,  the 

n  h 

total  variability  is  fi,    -  0.45/h  .  Variability  in  the  width  of  beams  is  quite  small  [41]. 

h  n 


Basic  variability  6,  in  the  placing  of  reinforcement  is  illustrated  in  Figure  A. 3. 
d 

As  with  member  thickness,   6,  is  seen  to  decrease  with  member  size.     Measurements  of  effective 

d 

depth  to  reinforcement  in  slabs  [-32,33]  were  taken  directly.     In  Drysdale's  study  [18]  of 
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Figure  A. 3  -  Variability  in  Effective  Depth  to  Reinforcement 
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columns,  however,  statistics  for  the  ratio  r  =  d'/d'   ,  the  ratio  of  the  actual  to  nominal 

n 

concrete  cover,  were  reported  instead  for  column  sizes  ranging  from  12  in  (305  mm)  to  42 

in  (1067  mm)  in  size.     In  order  to  extract  &,  from  the  statistics  of  r,  note  that  d  =  h  - 

d 

d'  r,  assuming  symmetrical  reinforcement  in  the  column,  and  thus 
n 

"           +  (d')  V]^/^  (A. 5) 

d         _n  n  r  

h    -  d'  7 
n  n 

In  order  to  use  this  relationship,   it  was  assumed  that  d'/h    was  1/8  or  d'  =  2  in  (51  mm) 

n    n  n 

and  that  6^  -  0.4/h^  as  previously  determined;  the  resulting  6^  are  plotted  in  Figure  A. 3. 

It  appears  that  6,  -  0.6/h  ,  which  is  somewhat  higher  than  for  member  dimensions.  Birkland 
d  n 

and  Westhoff  [7]  suggest  that  variations  in  placement  ranging  from  +  1"  to  -1  3/8"  (+25 

to  -35  mm)  may  be  common  in  practice;  using  Equation  3.2b,  this  leads  to  an  estimate  of  6^ 

which  is  even  slightly  higher.      The  study  by  Mirza  and  MacGregor  [41]  suggests  that  6,  - 

a 

0.594/d  in  slabs, 
n 

Additional  variability  is  induced  by  the  deviation  of  mean  effective  depth  from  d^. 

The  ratio  d/d    appears  to  be  consistently  biased,  varying  between  0.97  and  1.0.    While  for 
n 

slabs,  this  A,  -  0.05  [32,33]  it  appears  to  be  about  0.02  for  columns  [18].    This  implies 
d 

that  A,  -  0.32/h  ;  the  total  variability  is  then  fi,  =  0.68/h  . 
d  n  d  n 

The  data  discussed  above  have  been  derived  for  concrete  members  which  are  cast  in 
situ.     Fiorato's  recent  study  [25]  also  contains  data  on  variabilities  in  precast  members. 
These  data  indicate  that  the  variability  in  precast  member  geometry  is  less  than  one-half 
that  for  members  cast  in  situ. 


A. 2    Structural  Resistance 
A.  2.1  Flexure 

The  calculated  ultimate  flexural  capacity  of  a  properly  designed  underreinf orced  concrete 
beam  according  to  ultimate  strength  provisions  of  ACl  318-71  is, 

A  f 

M    =  A  f  d  (1  -  0.59  rf  r^)  (A. 6) 

u        s  y  bd  f 

c 
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The  mean  and  c.o.v.  may  be  calculated  as, 


  A  f 

M    ^  A  f    d  (1  -  0.59  —  (A. 7a) 

^  y  bd  f 

c 

=   [(1  -  1.18  q)2 
U  y 


+  (0.59  q)^  fi^    +  S]^]/(l  -  0.59  q)^  (A. 7b) 

c 


using  Equation  3.3  in  which  q  =  (A  f  /bd  f  )  =  p  f  /f  .  Contributions  of  variability  A 

syc              yc  I 

and  b  have  been  neglected,  as  they  and  the  coefficients  by  which  they  are  multiplied  in 
Equation  A. 7b  are  quite  small. 


Additional  uncertainties  arise  from  the  inability  of  Equation  A. 6  to  predict  the 
failure  load  precisely  when  all  its  parameters  are  known  (modeling  errors) .     This  contri- 
bution of  bias  and  variability  is  analyzed  using  data  provided  by  Corley  [13]  and  Sexsmith 
[55]  and  Equations  3.5.     The  values  B  =  1.11  and  o/B  =  0.07  are  based  on  laboratory  test 
data,  wherein  all  parameters  presumably  are  known  with  some  confidence,  and  thus  represent 
bias  and  variability  due  to  the  equation  itself. 

2 

The  total  variability  Q„  in  flexural  resistance  is  thus  evaluated  as  fi,,  =  [Q,,  + 

M  M 

(0,07)^]"'"''^.     With  Grade  60  reinforcement,  4000  psi  (2758  KN/m^)  concrete,  and  p  =  0.02, 

parameter  q  =  0.3.     With  these  values,        is  found  to  decrease  from  0.159  when  h  =  8  in 

(203  mm)  to  0.126  when  h  =  24  in  (610  mm).     In  Ref .  21         was  found  to  vary  slightly  with 

p,  but  this  variation  is  small  and  can  be  ignored. 


The  ratio  between  the  mean  in-situ  resistance  and  the  nominal  calculated  resistance 
(capacity  computed  using  nominal  material  strengths  and  dimensions)   is  also  required  for 
calibration.     Using  the  data  presented  in  Section  A.l,  M/M'   is  found  to  average  about 
1.12.     Any  reduction  in  concrete  strength  under  long-term  load  has  little  effect  on  flexural 
strength,  since  this  is  determined  primarily  by  A  f    and  d. 
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A.  2. 2  Shear 

The  analysis  of  uncertainties  in  shear  resistance  is  considerably  more  complex  than 
for  flexure.     The  truss  analogy  model  used  to  compute  the  capacity  is  approximate  in 
nature  because  the  mechanism  of  shear  in  beam  and  slabs  is  not  completely  understood. 
Recent  studies  [54,61]  have  shown  the  c.o.v.   in  shear  resistance  to  be  in  the  range  0.18 
to  0.24.     Because  of  the  implications  that  such  high  variabilities  have  for  structural 
safety  calculations,   the  analysis  of  uncertainties  in  shear  is  discussed  at  some  length  in 
the  following. 

According  to  the  truss  analogy  equation,  the  ultimate  shear  capacity  of  reinforced 
concrete  beams  is  given  by 

V  =  V    +  V    =  v    bd+         A    f  (A. 8a) 

u        c        s        c  s      V  y 

or 

V  =  v    +  p     f  (A. 8b) 
u        c        V  y 

in  which  v    is  the  nominal  shear  stress  at  which  inclined  cracking  occurs; 
c 

v    =  1.9  /P"  +  2500  p  V  ^  <  3.5  /f^  (A. 9a) 

c  c  M  —  c 

and  p     =  A  /bs,      A    and  s      are  area  and  spacing  of  the  stirrups. 

V  V  V 

A  great  number  of  tests  of  reinforced  concrete  beams  in  shear  conducted  during  the 
past  twenty  years  are  summarized  in  reports  of  the  joint  ACI-ASCE  Committee  426  [53,61]. 
The  variability  in  these  tests  is  illustrated  in  Figure  A. 4,  where  test  values  of  v^  are 
plotted  against  the  v^  calculated  from  Equation  A. 8b.  As  these  tests  were  conducted  under 
laboratory  conditions,  the  variability  may  be  interpreted  as  that  which  would  exist  when 
all  parameters  in  the  prediction  model  are  known.  Inspection  of  this  large  body  of  test 
data  reveals  that  Equation  A. 8b  is  conservatively  biased. 
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Figure  A. 4  -  Bias  and  Variability  in  Shear  Equation 
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A  portion  of  this  bias  and  variability  is  caused  by  Equation  A. 9a,  which  is  also 

empirical.     Table  5.20  of  ref.  53  summarizes  the  mean  and  c.o.v.  of  the  ratio  v  /v 

c,test  c,calc 

for  430  beams  without  web  reinforcement.  In  that  analysis,  data  with  varying  characteristics 
were  lumped  together.     In  an  effort  to  obtain  a  more  homogeneous  sample,  these  data  were 


first  separated  according  to  source  and  to  whether  v  /bd  /f  was  less  than  or  equal  to 

c  c 

3.5.     The  mean  and  standard  deviation  of  the  ratio  of  test  to  calculated  v    were  then 

c 

recomputed.     The  results  are  summarized  in  Table  A.l.     The  mean  and  variance  viewed  on 

this  basis    both  tend  to  be  lower  than  those  values  reported  in  Table  5.20  of  ref.  53. 

—  2  2 

Statistically  unbiased  estimates  of  this  mean  and  variance  are  B  =  1.105  and  o    =  (0.116)  . 
Inclusion  of  the  more  recent  data  of  Rajagopolan  [49]  and  Kani  [36]  does  not  change  these 
estimates  significantly. 

Additional  bias  and  uncertainty  arise  from  the  truss  analogy  model  itself,  including 

its  failure  to  consider  shear  friction  and  doweling,  and  the  assumption  that  the  stirrups 

have  yielded.     However,  it  is  difficult  to  sort  out  these  various  sources  of  bias  from 

available  data  since  the  actual  portions  of  shear  carried  by  the  concrete  and  steel  in  a 

beam  with  web  reinforcement  are  unknown.     Much  of  the  data  presented  in  Tables  6.1  and  6.2 

of  ref.  54  are  not  helpful  in  answering  these  questions,  because  in  many  cases  the  limiting 

values  of  8  /V'  or  10  /V  were  used  to  calculate  the  capacity  and  the  contributions  of  the 
c  c 

steel  and  concrete  cannot  be  correctly  apportioned.     Moreover,  it  was  apparent  that  the 
equation  bias  is  not  constant  for  all  beam  geometries,  but  tends  to  increase  as  the  parameter 

p  f  /v    increases.     Finally,  the  bias  tends  to  be  higher  for  restrained  beams  and  for  T- 

V  y  c 

beams. 

Data  on  beams  with  web  reinforcement  tested  in  shear  [46,54]  were  separated  according 
to  source  and  whether  v^  <  8  /f^  in  an  attempt  to  obtain  more  homogeneous  samples  than 
those  analyzed  in  Tables  6.1  and  6.2  of  ref.  54.     Beams  with  diagonal  stirrups,  those 
which  failed  by  splitting,  and  T-beams  were  omitted.     The  analysis  of  these  data  is  summarized 
in  Table  A. 2.     The  c.o.v.   in  the  ratio  of  test  to  calculated  strength  tends  to  be  lower 
than  the  values  0.2  and  0.22  presented  previously  [54]. 


The  simplest  approach  for  incorporating  bias  and  variability  arising  from  the  truss 

2 


analogy  equation  is  to  apply  an  overall  factor  B  to  V^.     The  quantities  and 
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Table  A.l  -  Analysis  of  Inclined  Cracking  Equation 


Source  n  3  o 


Moody,  et.  al  (S)* 

28 

1 

164 

0 

125 

u-i 

Moody,  et.  al  (R)* 

15 

1 

082 

0 

140 

ro 
V 

Bower  (S)* 

5 

1 

004 

0 

090 

Bower  (R)* 

28 

1 

004 

0 

090 

Rodriquez  (c)* 

13 

1 

001 

0 

109 

^" 

Thurston  (S)* 

58 

1 

241 

0 

129 

PCA  -  anson* 

5 

1 

196 

0 

063 

Rajagopolan  (S)[49] 

8: 

1 

144 

0 

208 

O 

> 

Kani  (S)  [36] 

13 

1 

107 

0 

141 

II  m 

Moody,  et.  al.*  (S) 

12 

0 

884 

0. 

094 

Moody,  et.  al.*  (R) 

43 

1 

052 

0 

112 

d 

R  =  restrained;  S 

=  simple  beam 

y 

> 

*Data  from  Chapter 

5,  ACI  Committee 

326  Report 

[54]. 

Table  A. 2  -  Analysis  of  Truss  Analogy  Equation 


Source  n                             3  o 

Clark  (S)*  37  1,325  0.229 

Bresler  (S)*  9  1.475  0.306 

r           Thurston  (S)*  19  1.165  0.121 

00           Placas  (S)   [46]  17  1.421  0.159 

V           Elstner  (R)*  4  1.549  0.122 

3         Rodriguez  (R)*  8  1.045  0.106 


I 

00 


It 


Clark  (S)*  9                          1.198  0.106 

Elstner  (R)*  23                           1.601  0.227 

Rodriguez  (r)*  12                          1.523  0.304 

*Data  from  Chapter  6,  ACI  Committee  326  Report [54] . 

R  =  Restrained  beam;  S  =  simple  beam 
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this  approach  would  be  determined  from  linear  statistical  analysis  of  Equation  A. 8a,  while 
2  2  2 

would  equal  V^o  ,  conforming  to  the  increase  in  scatter  observed  in  Figure  A. 4  as  v 

—  2 

increases.     Statistically  unbiased  estimates  for  B  and  o  ,  obtained  from  the  data  in 

Table  A. 2,  are  B  =  1,31,        =  (0.197)^  and  a/B  =  0.150  if  v    <  8  /T.     This  simple  model 

u  c 

can  include  terms  where  v    =8  /F^  (or  10  /F^)  in  the  equation  uncertainty  analysis,  as  we 

u  c  c 

are  looking  simply  at  overall  bias.     However,  when  P^f^  becomes  small,  approaching  the 

case  of  a  beam  without  web  reinforcement,   the  equation  bias  should  be  about  1.11  in  accordance 

with  the  previous  analysis  of  the  inclined  cracking  load  equation.     This  approach  thus  is 

incompatible  with  the  data  which  show  that  the  bias  increases  as  p  f  /v  increases. 

v  y  c 

To  eliminatii  this  incompatibility,  individual  bias  terms  may  be  applied  insteaa  to  tne 
portion  of  shear  carried  by  the  concrete  and  web  steel.  ,  The  mean  and  variance  would  be 


V  =  B,v    bd  +  B„  -4:-  A  f  (A. 10a) 

1  c  2     s      V  y 

2—2     2  —2     2  2 

c  s 

2  2  2  2 

If  a    =(v    +pf)a,a    may  still  be  interpreted  as  the  variance  in  the  ratio  of  test 
e  c        V  y 

to  calculated  strength.     So  that  Equation  A. 10  is  valid  when  p  f  /v    =  0,  B, ,  is  first  set 

V  y    c  1 

equal  to  1.11;  B^  =  1.38  is  then  determined  by  least  squares  fitting  of  that  data  (Table 

A.  2)  where  v    <  8  /F"  (or  10  /P") . 

u  c  c 


The  total  variability  in  resistance  to  shear  is  9,^  =  o^/V,  were  V  and        are  obtained 

from  Equations  A. 10  and  and          are  determined  from  linear  statistical  analysis  of 

c  s 

Equation  A. 8.     The  appearance  of  the  term  /f^  in  Equation  A. 9a  is  intended  to  reflect  the 

tensile  strength  f^  of  the  concrete.     Thus,  by  replacing  /f^  with  x  f^,  where  x  is  a 

constant,  it  is  found  that        ^9..    =  Q,^  .     This  result  also  can  easily  be  shown  to  obtain 

v        f  f 
c        t  c 

in  those  situations  where  v    =3.5  /f^,  or  at  very  small  reinforcement  ratios  where  it  has 

c  c' 

been  suggested  [61]  that  the  inclined  cracking  equation  should  be 


v    =  (0.8  +  120p)  /P"  <  2.3  /P"  (A. 9b) 

c  c  —  c 
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While        tends  to  be  slightly  higher  for  very  lightly  reinforced  webs,  for  practical 

purposes     it  may  be  treated  as  constant  at  Q„  =  0.22.     The  variation  in        and  caused 

V  d  b 

by  member  size  has  little  effect  in  Q.^,  because  of  the  magnitude  of  the  other  variabilities, 
particularly  that  due  to  the  equation. 

Calculations  using  Equations  A. 10  and  the  data  in  section  A.l  show  that  V/V'  ranges 
from  about  1.14  when  the  minimum  web  reinforcement  permitted  by  ACT  318-71  is  used  (A  = 
50  bs/f^)  up  to  about  1.30  for  very  heavily  reinforced  webs.     For  purposes  of  calibration 
to  existing  designs,  a  value  of  1.18  will  be  used,  recognizing  that  the  results  may  be 
slightly  conservative  for  heavily  reinforced  webs. 

A. 2. 3    Axial  Thrust  and  Bending 

The  analysis  of  uncertainty  in  the  resistance  of  reinforced  concrete  sections  subjected 
to  a  combination  of  axial  thrust  and  bending  presents  special  problems.     This  resistance, 
denoted  R  ,   is  described  by  an  interaction  diagram,  a  locus  of  thrust  and  moment  combinations 
(P,  M)  that  cause  failure.     General  expressions   [45]   for  resistance  derived  directly  from 
ultimate  strength  principles  are  unavailable  in  closed  form.     Therefore,  Equation  3.4 
cannot  be  applied  in  evaluating  the  mean  and  variance  in  resistance. 

Monte  Carlo  simulation  provides  a  useful  alternate  method  for  evaluating  the  statistics 
R  .     Values  of  f   ,   f  ,   etc.  are  sampled  numerically  from  probability  distributions  having 

u  ye 

means  and  c.o.v.  described  in  section  A.l  and  a  series  of  interaction  diagrams  is  generated. 

The  mean  and  c.o.v.  of  R    can  then  be  evaluated  as  a  function  of  the  eccentricity  ratio 

u 

e/h  =  M/Ph  or  at  particular  mean  thrust  or  moment  levels,  depending  on  the  character  of 
the  applied  load  [24,52].     If  thrust  and  moment  increase  in  proportion  to  failure,  the 
statistics  at  fixed  eccentricity  are  most  appropriate. 


This  analysis  has  been  described  in  some  detail  in  a  recent  publication  [22] .  Therein, 

2  2  1/2 

6„     and  A„     were  analyzed  separately;   according  to  Equation  3.1,  =   [&      +  '^n  ] 

R  R  R  R  R 

u  u  u  u  u 

f2„    may  also  be  obtained  directly  using  the  estimates  of  total  variability  (Q,  ,Q   ,  etc.) 
R  b  n 

u 

presented  in  Section  A.l  in  the  simulation.     Q,^     is  shown  as  a  fuction  of  p     and  e/h  in 

^u  § 

84 


Figure  A. 5  for  a  short  beam-column  constructed     of  4000-psi  concrete  and  Grade  60  rein- 
forcement with  average  quality  control.     Simulations  on  sections  of  different  size  yielded 

essentially  the  same  results.     When  p     is  small,  Q      is  elevated  at  small  e/h;  however,  it 

S  ^u 

decreases  as  the  thrust  component  lessens  and  the  compressive  strength  of  the  concrete 

becomes  less  significant  to  the  overall  resistance.     It  appears  that  il      may  be  taken  as 

R 

u 

sensibly  constant  provided  that        >  0.02  (22).     Thus,     except  for  lightly  reinforced 
columns,  the  increase  in  ((>'  with  e/h  permitted  in  ACI  Standard  318-71  must  be  justified  on 
the  basis  of  the  accompanying  additional  ductility  prior  to  failure  rather  than  any  significant 
decrease  in  fi^.     As  a  point  of  interest,  the  Monte-Carlo  simulation  should  yield  approximately 
the  same  results  as  Equation  A.  7b  as  e/h  ->■  «>.     With  p  =  0.03,  the  latter  equation  yields 

=  0.125,  which  may  be  seen  to  compare  quite  well  with  the  simulated  values  in  Figure  A. 5. 

u 

The  bias  and  variability  in  the  strength  prediction  are  assessed  from  analyzing 

—  2  2 

laboratory  tests  of  short  columns  [31].     These  data  show  that  B  =  0.99  and  a    =  (0.0603)  , 

making  the  contribution  of  equation  variability  to  Q    about  0.06.     That  this  is  close 

R 

to  the  flexural  equation  variability  is  explained  by  the  similarity  in  the  assumptions 

underlying  the  flexural  and  combined  thrust  and  bending  analyses  [39,45].     Using  the 

results  in  Figure  A. 5,  most  columns  would  fall  in  the  range  0.14  <  Q    <  0.18.     The  value 

R 

of  0.17  used  in  the  calibration  studies  corresponds  to  a  moderately  reinforced  column 
constructed  with  average  quality  control. 

The  ratio  between  the  mean  in  situ  and  nominally  calculated  resistance,  calculated 
using  methods  similar  to  those  for  flexure  and  shear,  was  found  to  depend  on  e/h.     A  best 
overall  value  was  0.965.     This  is  lower  than  for  flexure,  since  the  difference  between  in 
situ  and  mean  concrete  discussed  earlier  becomes  more  important  for  sections  which  must 
resist  significant  axial  loads. 

The  data  used  above  to  analyze  the  equation  bias  and  variability  were  obtained  from 
tests  on  tied  columns.     The  results  should  not  differ  significantly  for  columns  with 
spiral  reinforcement,  since  the  primary  advantage  in  using  spiral  reinforcement  is  in  the 
greatly  improved  column  ductility  that  it  affords  [38,39]. 
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Figure  A. 5  -  Coefficient  of  Variation  in  Beam  Column  Resistance 

Due  to  Coefficients  of  Variation  in  Individual  Variables 

(f  =  60000  psi,  f  =  4000  psi) 
y  c 
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Slenderness  effects  in  columns  are  accounted  for  in  existing  designs  using  the  moment 

magnifier  method  which  amplifies  the  design  moment  on  the  column.     This  may  affect  the 

bias  and  variability  in  R.     In  an  extensive  comparison  of  test  values  to  predictions  for  a 

kJi/r  <_  100  it  was  found  [38]  that  E(Test/Calc)  =  1.13,  with  a  c.o.v.  of  0.169.     These  data 

are  stratified  according  to  k!L/r.     Analysis  of  those  data  where  kH/r  =  46  -  57,  for  example, 

leads  to  a  =  0.152.     Removing  that  component  of  variability  inherent  to  the  short  column 

analysis  (0.06),  the  c.o.v.  attributable  to  the  moment  magnifier  method  is  about  0.13  for 

columns  where  k£/r      70.     Studies  have  shown  [38]  that  k£/r  usually  is  less  than  70  in 

reinforced  concrete  columns;   in  fact,   in  most  braced  frame  structures,  kl/r  <  40.  Thus, 

for  slender  columns,        would  increase  from  0.17  to  about  0.21,     In  terms  of  safety, 

R 

however,  this  is  partially  offset  by  the  effect  of  the  more  conservative  equation  bias  (B 
-  1.13)  on  R/R' . 

A. 3    -  Structural  Loads  and  Load  Effects 

Variabilities  in  load  effects  which  arise  from  gravity  and  wind  are  emphasized  herein 
since  for  most  structures  they  are  the  loads  often  checked  in  preliminary  designing.  The 
basic  approach  taken  is  to  relate  the  load  A  to  its  load  effect  S^  using  the  influence 
coefficient  c. : 


assuming  the  transformation  from  load  to  load  effect  is  a  linear  one.     The  total  variability 

is  defined  by 


in  which  uncertainties  in  c^  arise  from  the  analysis  transforming  the  load  to  a  load 
effect  and        would  include  basic  variability  in  A,  as  well  as  uncertainty  in  the  load 
model  itself.     It  may  be  noted  that  the  term  Q      in  Equation  A. lib  serves  the  same  purpose 


A' 


(A. 11a) 


(A. lib) 


as  the  "analysis"  term  proposed  by  Cornell  [15]    and  used  by  Galambos[28] . 
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A. 3.1    Gravity  loads 

The  dead  or  permanent  load  results  from  the  weight  of  elements  comprising  the  structure 
and  also  includes  permanent  equipment,  partitions  and  installations.     This  weight  is  quite 
predictable  when  the  geometry  of  the  structure  is  specified  and  equipment  specifications 
are  known.     The  ratio  of  mean  to  nominal  value  is  assumed  to  be  unity.     Previous  studies 
have  employed  values  of        ranging  from  0.06  to  0.20  [21,24,28,37],  most  being  between 
0.06  and  0.13;  herein,  a  value  of  0.10  is  used. 

The  specification  of  live  loads  for  structural  design  consists  of  a  basic  uniformly 
distributed  unit  load  (load  per  unit  area) ,  which  depends  on  a  particular  occupancy  type, 
and  a  reduction  factor  by  which  the  basic  load  may  be  multiplied  to  account  for  the  decrease 
in  unit  load  observed  as  the  loaded  or  tributary  area  increases.     The  load  reduction 
factor  currently  recommended  by  ANSI  A58. 1-1972  and  other  standards  can  be  traced  back 
thirty  years  where  it  was  determined,  somewhat  subjectively,  from  load  survey  data  gathered 
in  two  federal  office  buildings. 

In  recent  years,  a  significant  amount  of  load  survey  data  have  become  available  for 
office  buildings  in  the  United  Kingdom  (UK)    [42]  and  the  United  States  (US),  the  latter 
data  being  gathered  in  a  recent  NBS  load  survey  [16] .     Concurrently,  analytical  procedures 
to  model  live  loads  on  structures  probabilistically  [40,44]  have  been  developed.  The 
methods  enable  an  equivalent  uniformly  distributed  load  (EUDL)  to  be  determined  which  will 
produce  (statistically)  the  same  load  effect  on  a  structural  member  (beam  moment,  column 
thrust,  etc.)  as  the  actual  random  set  of  loads.     This  methodology  has  been  applied  successfu 
to  the  UK  load  survey  data  [40],  and  to  the  US  data  [23],  as  summarized  here. 

On  the  basis  of  a  statistical  analysis  of  NBS  data  from  all  rooms  which  were  randomly 

2 

selected,  the  total  mean  unit  live  load  was  found  to  be  11.6  psf  (555  N/m  );  this  is  also 
equal  to  the  mean  EUDL.     In  comparison,   the  mean  unit  load  obtained  using  the  UK  data  [40] 
was  11.8  psf  (565  N/m^) .     The  variance  in  L  was  found  to  be  Var  [L]  =  26.2  +  14300/A,  in 
which  A  is  the  influence  area. 

The  EUDL  L  obtained  from  load  survey  data  L  denotes  the  point-in-time  sustained  live 

load.     The  probability  distribution  of  the  maximum  sustained  load  L^  seen  during  the 
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lifetime  T  of  the  building  is  determined  by  assuming  that  occupancy  changes  occur  as  a 

Poisson  random  process  with  rate  v  ;  the  distribution  of  L    is  then  given  by, 

Ij  s 

(x)  =  exp  [-Vj^Td  -  F^(x))]  (A. 12) 

s 

in  which        is  the  distribution  function  for  L.     NBS  survey  data  show  an  average  occupancy 

duration  of  8  years;  with  T  about  64  years,  v^^T  =  8.     Loads  may  also  occur  from  infrequent 

clustering  of  people  above  and  beyond  normal  personnel  loads  and  other  activities  such  as 

remodeling  (extraordinary  loads) .     For  the  purposes  of  structural  design,  the  mean  lifetime 

maximum  total  EUDL        is  of  interest,  as  it  is  this  EUDL  which  is  analogous  to  the  design 

live  load  specified  in  building  codes  and  standards.     The  upper  fractiles  (0.9  -  0.99)  of 

may  be  estimated  as  the  corresponding  upper  fractiles  of  the  maximum  sustained  load 

plus  the  mean  of  the  extraordinary  load  for  the  period  over  which        acts   [40] .  Approximate 

formulas  for  the  mean  and  variance  of  L    or  L    may  be  derived  by  fitting  a  Type  I  Extreme 

t  s 

Value  distribution  to  the  upper  load  fractiles  and  calculating  the  mean  and  variance  of 
the  fitted  Type  I  distribution. 


A  comparison  of  expressions  for  E[L^]  and  E[L^]  with  the  live  loads  currently  permitted 
by  ANSI  A58. 1-1972  is  presented  in  Figure  A. 6.     As  the  loads  are  presented  as  functions  of 
influence  area  rather  than  tributary  area,  separate  curves  are  required  for  ANSI  floor  and 
column  loads.     The  agreement  between  values  of  E[L^]  based  on  the  UK  and  NBS  data  is  very 
close.     It  may  also  be  observed  that  the  maximum  sustained  load  is  much  less  than  the  ANSI 
load  where  areas  are  small,  but  is  almost  identical  to  it  at  larger  areas. 


For  the  maximum  total  EUDL  L^  the  analysis  described  above  gives,  approximately, 

E[L^]  =  18.7  +  520//A  (A. 13a) 

Var[L^]  =  14.2  +  18900/A  (A. 13b) 

The  analysis  using  the  UK  data  resulted  in  E[L^]  =  14.9  +  763//A  and  Var  [L^]  =  11.3  + 
15000/A.     As  seen  in  Figure  A. 6  there  is  a  significant  difference  between  the  two  mean 
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loads  at  small  influence  areas.     This  difference  is  attributed  to  differences  in  the 
parameters  used  in  the  extraordinary  load  model.     This  is  somewhat  disturbing,  as  unfortunately 
there  is  little  objective  basis  for  selecting  the  extraordinary  load  parameters.  Additional 
research  is  needed  so  that  the  parameters  of  the  extraordinary  live  load  model  can  be 
estimated  with  confidence. 


At  large  areas,  E[L^]  exceeds  the  ANSI  live  load  regardless  of  which  extraordinary 

load  model  is  chosen.     It  seems  clear  that  the  maximum  load  reduction  of  60  percent  allowed 

by  ANSI  is  too  large;  according  to  this  analysis,  the  maximum  load  reduction  should  be 

something  less  than  50  percent.     Moreover,   the  current  load  reduction  rate  of  0.08  percent 

2  2 

per  square  foot  of  area  over  150  ft     (13.94  m  )  permits  too  rapid  a  reduction  in  load  for 
beam  loads,  although  it  appears  to  be  about  right  for  column  loads. 


The  total  variability  fi^  in  maximum  live  load  effect  is  obtained  by  augmenting  the 
data-based  variability,  obtained  as    Var [L^ ] /E[L^ ]  from  Equations  A. 13,  with  uncertainty 
measures  attributable  to  the  load  model    and  the  structural  analysis  relating  load  and 
load  effect.     These  are  taken  as  0.20  and  0.05,  respectively  [15,21,28].     Q.^  thus  ranges 
from  about  0.28  for  A  =  200  ft^  to  about  0.26  for  A  =  2000  ft^,  and  for  purposes  of  this 


study  will  be  taken  as  constant  at  0.26. 


Statistics  for  the  instantaneous  live  load  are  required  when  combining  loads  from 

2  2 

several  sources  [63].     As  above,  this  uncertainty  is  obtained  by  adding  (0.20)     +  (0.05) 

1/2 

to  the  data-based  variability  which  is  given  by  [(26.2  +  14300/A) /II. 6]       .     Q  ranges 

^I 

from  0.71  when  A  =  400  ft^  (37. 2m)^  to  0.51  when  A  =  4000  ft^  (372m^) ,  a  variation  that 
must  be  considered  in  the  calibration. 


A. 3. 3    Wind  Loads 

The  wind  pressure  on  a  structure  is  given  by  an  expression  of  the  type 


W  =  GF[l/2p  C  (Tn)^"] 
p  iU 
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using  the  power  law  for  wind  velocity  profile  with  height  z,  in  which  p,  C    and  V  are  mass 

P 

density  of  air,  pressure  coefficient,  and  basic  wind  speed  at  a  height  of  10  m  (30  ft) 
respectively,  GF  is  the  gust  factor  which  reflects  the  dynamic  characteristics  of  the 
structure  and  which  enables  the  wind  forces  to  be  used  in  a  static  structural  analysis 
[17,64].     The  wind  force  W  may  be  written  as, 

W  =  W(Cp,V,K,a,f^,6)  (A. 15) 

in  which  k  and  a  are  exposure  terms  while  f^  and  B  are  fundamental  natural  frequency  and 
damping,  respectively.     Inspection  of  Equation  A. 14  reveals  that  the  major  component  of 
variability  in  W  arises  from  uncertainty  in  the  wind  speed,  V. 


When  gravity  and  wind  loads  are  combined,  account  must  be  taken  of  the  reduced  probability 
that  extreme  live  and  wind  loads  occur  simultaneously.     For  reliability  analyses,  the  load 
combination  of  interest  is, 

Z  =  max[L(t)  +  W(t)]  (A. 16) 

T 

where  T  is  the  life  of  the  structure.     The  duration  of  structurally  significant  winds  is 

almost  instantaneous  compared  to  that  for  live  loads.     Under  these  conditions,  McGuire 

and  Cornell  [40]   found  by  numerical  simulation  that  the  upper  percentiles  of  the  maximum  of  the 

sum  of  two  random  processes,  one  with  duration  and  one  without,  was  well  approximated  by 

the  same  percentiles  of  one  of  the  following: 


{L      +  W 
s  . 


(A. 17a) 
(A. 17b) 


where        and  L    are  as  defined  in  A. 3. 2,  W    is  the  maximum  of  W(t)  during  the  period  in 
Is  s 

which  L    acts  (here  about  8  yr)  and  W    is  the  lifetime  maximum  of  W(t) .     Analyses  by  Ferry 
s  T 

Borges  [24]  and  Wen  [65]  lead  to  similar  conclusions.     Accordingly,  the  analysis  of  reli- 
ability under  gravity  plus  wind  loads  requires  that  these  two  load  cases  both  be  con- 
sidered.    However,  studies  using  the  more  exact  procedure  suggested  by  Wen  [65]  show  that 

the  mean  of  Z  is  well  approximated  in  many  cases  by  the  mean  of  Equation  A. 17b. 

92 


Recent  statistical  analysis  of  wind  speed  data  conducted  by  the  National  Bureau  of 
Standards  [56]  shows  that  for  extratropical  winds  the  Type  I  Extreme  Value  distribution 
best  models  the  annual  maximum  wind  speeds.     With  this  distribution,  the  mean  and  c.o.v. 
of  the  maximum  V  of  n  identically  distributed  and  statistically  independent  variables  V 
are  related  by 


V  =  V     (1  -        6„    £n  -)  (A.  18a) 

o  IT    V  n 

o 

■p-  =  7^  -  —  £n  -  (A. 18b) 

<S^      6^       TT  n 

o 

Using  data  reported  by  Simiu  [56]  on  37  year  records  of  wind  speeds  excluding  tropical 

sites,   it  was  found  that  0.07  <  6^    <  0,24.     The  corresponding  range  in  c.o.v.  for  a 

o 

period  of  64  years  (equivalent,  on  the  average,  to  8  tenancies;  see  previous  section)  is 
0.06  <  6^  <  0.14,  with  0.11  as  a  representative  value;   for  8  years,  0.06  <  6^  <  0.17. 


The  mean  lifetime  wind  speed  may  be  related  to  the  wind  speeds  currently  specified  in 
codes,  e.g.  ANSI  A58. 1-1972,  which  are  based  on  mean  recurrence  intervals.  Corresponding 
to  a  mean  recurrence  interval  N,  the  wind  speed        is  (for  Type  I  distribution  for  V^)  is 

v^  =  V     [1  -  (S„     (0.45  +  0.7797  Sin  Sin  —  -)]         (A. 19) 

*        °  o  1  -  1/N 

Using  this  v^  and  V  defined  by  Equation  A. 18a,  it  may  be  seen  that  for  most  6^    and  N,  n 
of  interest,  V/v^  is  typically  1.05  -  1.07;   in  particular,  if  n  =  N,  V/v^  is  very  close  to 
1.0. 


Additional  uncertainty  in  the  wind  force  W  arises  from  variability  in  the  pressure 
coefficient  C  .     For  low-rise  structures,  fi„     falls  in  the  range  of  0.10  -  0.20  [20]  and 
may  be  even  higher  for  buildings  with  unusual  shapes.     For  more  conventional  buildings, 
this  uncertainty  is  somewhat  less,  and  a  value  of  0.15  is  used  in  this  study. 


The  total  variability  in  W  is  evaluated  as 
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„2  _  1      ^     9W    2  -2  2 

-  =2     "  ^3^.) 

W  1  1 


(A. 20) 


where        are  the  parameters  in  Equation  A. 14.     These  derivatives  have  been  evaluated  in 
Appendix  D  of  ref.   21  using  Davenport's  [17]  gust  factor.     The  gust  factor  analysis  presented 
by  Simiu  [57]  was  also  used  as  a  check.     For  example,  corresponding  to  the  wind  speed 
term, 


1     9W  -      „  , 

=    3V  ^  =  2  +  a^ 

w 

1     3 GF   

in  which  a    =  —  — —  V.     Using  Simiu 's  analysis,  a„  is  of  the  order  0.1-0.3.     In  contrast, 

v      —      9V  ^  V 

ref.  21  showed  that  0.1  <        1  0.4  for  buildings  ranging  from  100  -  600  ft  (30.5  -  183  m) 
in  height.     Thus,   the  contribution  of        to  fi^  may  be  obtained  by  multiplying  Q.^  by  roughly 
2.3. 


The  terms  k  and  a  are  related  to  the  building  exposure.     While  Q,    =0.20  and  - 

K  a 

0.10  [64],   the  coefficients  associated  with  them  in  Equation  A. 18  are  roughly  0.32  and 
0.4,  respectively  [21].     The  contribution  of  exposure  to  overall  wind  variability  thus  is 
small.     Finally,  the  natural  frequency  and  structural  damping  contribute  uncertainty  but 
because  the  derivative  terms  are  so  small  their  effect  can  be  neglected  in  a  practical 
sense.     The  total  variability  in  lifetime  maximum  wind  load  is  thus. 


n      =  [(2.3  X  0.11)^  +  (1  x  0.15)^  +  [(.32  x  .2)^  +  (.4  x  .1)^]  +  (0 . 05)  ^  ] ''■''^ 
T 

velocity  Pressure  Exposure  Modeling 

Coef f . 

a,    =  0.31 

T 

Of  this,  0.29  is  due  simply  to  the  velocity  and  pressure  coefficient  terms  which  are 

basically  structure- independent .     Similarly,  =  0.35. 

W 


The  ratio  of  riominal  to  mean  wind  loads,  required  for  the  calibration,  depends  on  the 
ratio  (v^/V),  the  ratio  between  the  ANSI  pressure  coefficient  and  the  mean  pressure  coef- 
ficient, and  the  orientation  of  the  building  to  the  direction  of  the  extreme  wind.  As 
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previously  noted,  v^/V  is  about  0.93  -  0.95.     The  pressure  coefficients  used  in  ANSI 

A58. 1-1972  tend  to  be  conservative  for  the  leeward  side  of  a  buildine;  rather  than  C  = 

-0.5,  the  mean  would  be  more  like  -0.3.     This  would  make  the  total  mean  drag  coefficient 

C'    ^  1.1  rather  than  1.3.     Data  contained  in  BSI  CP3:  Chapter  V:  Part  2:  1972   [11]   tend  to 
P 

support  this  position.     Finally,  there  is  a  reduced  probability  that  the  wind  will  act  on 
the  structure  in  its  most  unfavorable  direction.     The  effect  on  W  may  be  estimated  by 
assuming  the  orientation  of  the  force  vector  to  be  uniformly  distributed  between  -rr/A  and 
Tr/4  with  respect  to  the  axis  of  the  building,  leading  to  a  mean  force  which  is  from  85  - 
90  percent  of  the  most  unfavorable.     Combining  all  these  factors  leads  to 


W'/W^=  (0.94)^  X  (1.3/1.1)  X  (1/0.85)  =  1.2 

Similarly,  W / W    ^  1.75. 
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